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DEVELOPMENT  AND  HELD  TEST  OF  THE  CENTRAL  ENERGY  PLANT 
ADAPTIVE  CONTROL  SYSTEM  (CEPACS) 


1  INTRODUCTION 


Background 

U.S.  Army  installations  operate  and  maintain  a  large  number  of  fossil  fuel-powered  central  energy 
plants  (CEPs).  A  general  problem  at  Army  CEPs  is  that  outdated  plant  control  systems  have  become  both 
economically  and  operationally  inefficient  when  compared  to  modem  systems.  A  reliable  automatic 
control  system  is  an  essential  element  in  maintaining  safe,  efficient  CEP  operation.  The  difficulty  of 
acquiring  spare  parts  for  these  old  systems,  and  the  recent  price  reduction  in  modem  control  hardware 
and  software  have  made  it  more  practical  to  replace  the  obsolete  or  primitive  control  systems  used  in 
many  Army  facilities  with  controls  that  incorporate  state-of-the-art  electronics  and  advanced  boiler-control 
techrologies  to  improve  safety,  reliability,  and  efficiency. 

The  U.S.  Army  Construction  Engineering  Research  Laboratories  (USACERL)  and  the  University 
of  Illinois  Mechanical  Engineering  D^artment  are  developing  and  testing  a  prototype  advanced, 
self-tuning  Central  Energy  Plant  Adaptive  Control  System  (CEPACS)  for  retrofit  to  old  systems  or 
installation  in  new  boiler  applications.  This  prototype  system  promises  to  reduce  the  need  for  skilled 
operators  since  it  C2q)tures  plant  operating  experience  and  retunes  itself  for  optimal  performance.  This 
application  uses  a  control  strategy  based  on  the  Generalized  Predictive  Control  (GPC)  algorithm  (Clarke, 
Mohtadi,  and  Tuffs  1987),  which  has  been  widely  applied  in  a  number  of  process  industries  in  the  United 
Kingdom. 

Traditional  boiler  control  methods  base  their  control  schemes  on  the  general  assumption  that  plant 
dynamics  are  static  through  time.  The  GPC  algorithm  assumes  that  plant  dynamics  are  in  flux;  it  updates 
the  plant  model  parameters  after  each  sampling  period,  continuously  modifying  the  control  strategy  as 
needed  to  bring  the  plant  to  the  desired  set  of  conditions  in  the  minimum  amount  of  time.  The  model 
automatically  adjusts  itself  to  changes  in  equipment  configuration  or  behavior  as  well  as  to  ambient 
conditions.  Applying  the  GPC  to  boiler  control  is  a  new  concept  that  needs  to  be  field-tested  as  a  first 
stq)  toward  implementing  of  the  CEPACS  system.  This  field  test  of  CEPACS  was  designed  to  determine 
the  effectiveness  of  the  self-tuning  control  system  in  a  CEP  environment. 


Objectives 

The  objectives  of  this  project  were  to  develop  and  field  test  a  prototype  Central  Energy  Plant 
Ad^tive  Control  System  to  be  installed  in  Army  facilities  to  inq>rove  plant  safety,  reliability,  and 
efficimcy. 


Approach 

Lito-ature  on  boiler  control  basics  was  reviewed  and  analyzed.  The  prototype  self-tuning  adiq)tive 
control  system  was  developed,  including  the  control  algorithm  and  the  required  hardware  and  software. 
A  fidd  test  was  conduaed  on  a  gas-  and  oil-fired  boiler  at  the  University  of  Illinois  Abbott  Power  Plant 
in  July  1991 .  Test  results  were  analyzed,  and  further  system  enhancements  such  as  robustification,  faults 
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tolwance,  and  diagnostics  were  explored.  Such  essential  processes  as  the  automatic  startup/shutdown 
procedures,  furnace  draft  control,  and  on-line  calculation  of  the  boiler  combustion  efficiency  were 
inq)lemented  and  documented. 


Scope 

CEP  ACS  is  generally  applicable  to  gas-  or  oil-fired  boilers.  Solid  fuel  fired  units,  however,  require 
nKxlifications  of  the  control  strategies,  mainly  due  to  the  difficulties  in  accurately  measuring  and 
controlling  the  flow  rate  of  the  solid  fuel. 


Mode  of  Technology  Transfer 

It  is  recommended  that  the  information  presented  in  this  r^rt  be  incorporated  into  Design  Guide 
(DG)  and  an  Engineering  Technical  Note  (ETN). 
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2  BOILER  CONTROL  BASICS  AND  BOILER  MODELS 


Depending  on  system  complexity,  CEPs  use  any  of  several  different  types  of  control:  (1)  on/off 
control  (^so  called  two-position  control);  (2)  proportional  control  (also  called  gain  control  or  throttling 
control);  (3)  integral  control  (also  called  reset  control);  (4)  derivative  control  (also  called  rate  control);  and 
(S)  proportional-integral-derivative  control  (also  called  PID,  or  three-mode  control). 

The  use  of  conventional  proportional,  integral,  and  derivative  controllers  may  not  result  in  the 
highest  possible  operating  efficiency  because  of  the  large  variations  in  plant  inputs  and  disturbances 
affecting  the  process.  In  boilers,  these  variations  typically  arise  from  changes  in  fuel  calorific  value, 
warmup  effects  in  the  burner,  fouling  of  the  boiler  tubes,  or  deterioration  of  mechanical  links  and  profiles 
involv^  in  the  operation.  Climatic  changes  in  ambient  temperature,  pressure,  and  humidity  can  also  affect 
the  control  adversely.  For  the  control  system  to  deal  with  such  changes  effectively,  the  controller  must 
be  tuned  in  real  time,  preferably  without  operator  intervention.  This  can  be  done  by  implementing  a  loop 
between  the  control  signal  and  the  output  of  the  process  to  create  a  relationship  between  the  two.  This 
implies  the  generation  of  a  process  model  with  parameters  that  are  estimated  recursively.  These 
parameters  are  then  used  in  a  design  control  algorithm  to  compute  optimum  gains  for  the  controller.  The 
ability  to  estimate  model  parameters  in  real  time  for  a  range  of  loa^  and  inputs  results  in  online  tuning 
of  the  controller  gains. 

Since  self-tuning  is  faster  than  manual  tuning,  the  commissioning  time  for  installation  can  be 
decreased  by  as  much  as  50  percent.  Also,  self-tuning  performs  its  processes  systematically,  even  for  the 
simplest  control  loops.  Difficulties  have  been  encountered,  however,  with  processes  with  very  rapid 
parameter  variations,  or  very  strong  nonlinearities.  The  controller  cannot  be  applied  to  processes  that  do 
not  tolerate  the  process  se^int  deviation  required  in  the  identification  phase.  Difficulties  have  also  been 
found  under  operating  conditions  where  the  measured  value  is  suddenly  disconnected.  The  remedy  in  such 
cases  is  to  stop  the  parameter  updating  and  iwitch  the  controller  to  operate  in  the  conventional  Pn>  mode. 

PID  controllers  implemented  in  pneumatic,  electrcmic,  or  microcomputer  software  can  be  somewhat 
difficult  to  set  up  and  may  not  produce  stable  results  over  a  wide  range  of  operating  conditions.  Adaptive 
controllers  are  aimed  at  overcoming  these  deficiencies,  which  originated  from  dealing  with  control 
problems  related  to  high  performance  aircraft  and  rockets.  A  review  of  recent  adaptive  control  algorithms 
indicated  that  the  GPC  algorithm  was  the  algorithm  most  able  to  meet  the  control  requirements  of  CEPS. 
Researchers  developed  a  multivariable  mathematical  model  for  an  industrial  boiler  and  a  real  time 
simulator  for  a  steam  generation  system  to  evaluate  various  control  algorithms  and  operator  training.  This 
model  was  tailored  to  fit  an  actual  boiler,  and  the  tailored  model  was  then  used  to  redesign  the  original 
GPC-based  control  system.  The  system  also  contained  modifications  to  the  safety  interlocks  and  feed¬ 
forward  control,  and  used  the  identifier  information  extensively. 

In  enhancing  the  mathematical  model  of  the  boiler,  researchers  examined  many  analytical  models 
that  had  been  derived  to  predict  the  behavior  of  heating  plant  boilers.  Some  of  the  relevant  models 
recently  developed  were  compared  and  analyzed  to  determine  their  suitability  for  application  (Astrbm  and 
Bell  1988;  Bell  and  Astrbm  1987;  McDonald,  Kwatney,  and  Spare  1971;  and  Chawdry  and  Hogg  1989). 

Astrbm  and  Bell  (1987)  compared  several  previously  developed  models  (Morton  and  Price  1977; 
AstrOm  and  Eklund  1972,  1975;  Bell  and  AstrOm  1987),  ranging  from  simple  linear  to  nonlinear  models 
of  varying  degrees.  The  study  analyzed  the  most  effective  part  of  each  previous  model  and  combined 
them  into  a  single  nonlinear  model  that  used  a  combination  of  first  principles  and  field  data  to  determine 
model  parameters.  The  plant  studied  was  a  160  MW  oil-fired  natural  circulation  drum  unit  operated  by 
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Sydkraft  AB  in  Malmo,  Sweden.  The  plant  was  rated  at  1,108,800  IbAv  steam  at  1989  psig.*  The  inputs 
to  the  model  were  fuel  flow,  feedwater  flow,  and  control  valve  position.  The  plant  outputs  were  drum 
pressure,  drum  water  level,  steam  flow  rate,  and  electrical  output. 

This  nonlinear  model  structure  is  very  similar  to  the  required  form  of  the  research  system.  The 
electrical  output  subsystem  was  uncoupled  from  die  rest  of  the  equations,  so  separating  it  from  the 
remaining  part  of  the  system  was  simple.  The  oxygen  level  subsystem  could  be  added  easily  to  the 
equations.  Steam  flow  was  (tesignated  as  an  output  driven  by  state  equations  and  an  imaginary  value 
representing  the  net  sum  of  the  steam  use,  which  could  be  manipulated  into  the  test  boiler  configuration. 

To  suit  the  research  system,  several  modifications  and  additions  were  made  to  the  first  AstrOm  and 
Bell  model  (1987).  An  equation  governing  the  oxygen  level  loop  was  developed  and  added  to  the  system. 
The  equation  considers  a  stoichiometric  reaction  of  fuel  and  air  with  a  first-order  lag.  Additional  regime 
is  also  needed  during  the  startup  and  shutdown  periods,  where  the  boiler  passes  through  several  different 
regimes.  By  altering  the  constants  in  the  model  to  cover  a  wider  range  of  conditions  while  still  meeting 
the  steady  state  operation  model  at  the  operating  point,  a  general  model  to  go  from  the  start  to  the  steady 
state  can  be  obtained.  It  is  also  possible  to  account  for  the  fact  that  the  drum  is  circular.  For  normal 
operation,  the  behavior  of  the  water  level  is  close  to  linear,  but  if  the  water  level  approaches  the  extremes 
of  the  drum,  the  nonlinearities  become  significant.  A  solution  to  the  nonlinear  equation  is  very  difficult 
to  find,  but  by  fitting  a  polynomial  to  the  equation,  a  relationship  between  the  change  in  height  and  the 
change  in  volume  can  be  derived.  One  fin^  change  in  the  model  was  to  add  actuator  dynamics  to  the 
input  signals.  These  dynamics  can  be  modeled  by  a  pure  time  delay  and  a  computation  lag.  The  boiler 
simulation  models  in  state  space  structure  are: 
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t 


=Cu*l  C.iU.X,  +C,jU,  c„u, 


Cj.Uj  CjjU, 

Xj  =  -  CjjXj 

U,+Uj 


^^3  =  Cj.X,  CjjU, 


yi  =  ^51*1 


[Eq  1] 


(c^,Xj  +  l)(x  +C<j)  X,  X 

~  C45U2-C46)(“4C47*i) 

XjX, 


Where 

X,  is  the  pressure 
Xj  is  the  oxygen  level 
X3  is  the  fluid  density 
u,  is  the  fuel  flow 
Uj  is  the  air  flow 
U3  is  the  feedwater  flow 
U4  is  the  steam  flow. 


*1  lb  >  0.453  kg;  1  psi  «  6.89  kPa. 
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The  outputs  y„  yj,  and  are  drum  pressure,  oxygen  level,  and  water  level,  respectively.  All  the  c^’s  are 
the  proper  scaling  constants. 

For  control  purpose,  a  linear  boiler  model  is  needed,  and  can  be  reasonably  obtained  by  performing 
the  linearizatimi  of  the  nonlinear  model  around  the  boiler  operating  region.  The  boiler  model  in  the  linear 
state  space  form  is: 


*3  =  B33U3-^B34“4 

yi  =  c„x, 

y2  “  ^22^ 

y3  =  Cj.x,  +C33X,  4.D33U3  +  D3,u, 


[Eq2] 


Where  A,j’s,  By’s,  C,j’s  and  D^’s  are  the  constants  determined  by  the  experiments. 

The  test  boiler  was  a  dual-fuel  (oil  and  gas)  Bred  unit  enable  of  producing  175,000  Ib/hr  of  steam 
at  a  pressure  of  325  psig.  The  steam  from  this  no.  2  boiler  and  the  other  gas  boilers  flows  into  a  common 
hea(ter  with  the  steam  throttled  by  a  turbine  and  a  Republic  valve  in  normal  operatiem.  Researchers  were 
interested  in  controlling  the  header  pressure,  the  level  of  the  water  in  the  drum,  and  dte  oxygen  level  in 
the  flue  gas.  These  outputs  should  be  maintained  despite  variations  in  steam  demand,  fluctuations  in  the 
heating  value  of  the  fuel,  or  other  disturbances  to  the  system.  To  meet  these  control  objectives,  the  cmitrol 
system  presently  installed  is  capable  of  actuating  gas  flow,  oil  flow,  air  flow,  and  feedwater  flow  while 
sensing  steam  pressure,  steam  temperature,  steam  flow,  fuel  flow,  drum  water  level,  feedwater  flow,  air 
flow,  and  flue  gas  oxygen  level.  To  improve  the  boiler  model  to  match  the  actual  plant,  operating  data 
was  collected  from  the  plant  This  required  the  interface  of  the  developed  software  with  plant  equipment 
the  actual  data  collectirm  process,  and  data  manipulation.  It  further  required  altering  die  nonlinear 
equations  to  be  consistoit  with  these  results.  The  validity  of  the  final  model  was  verified  through 
simulations  (Appendix  A). 

Based  on  physical  principles  and  the  analysis  of  real  test  data  obtained  in  July  1991  at  the  Abbott 
boiler,  the  recommended  model  structures  for  each  loop  (the  model  order  and  corresponding  time  delay) 
are: 

1 .  Model  for  drum  pressure  loop: 

A(q  -')y (k)  =  q  -"B(q  -')u(k)  tEq  3] 

where 

A(q')  =  1  +  a,q'  + 

B(q')  =  bo  +  b,q'  +  bjq-2 
y(k)  is  the  drum  pressure 
u(k)  is  the  fuel  flow  command 
k  is  the  sample  number 
delay  n  =  10  seconds. 
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The  initial  parameters  used  for  this  loop  are: 


a,  =  -0.8959 
aj= -0.043 
bo  =0.2778 

b,  =  0.1503 
bj=  0.01523 

2.  Model  for  excess  oxygen  loop: 

A(q-')y(k)  =q-"‘B,(q-‘)u,(k)-^q-"'Bj(q-')u,(k) 


where 

A(q'')  =  1  +  a,q''  +  ajq'^  +  ajq  ^ 

B,(q')  =  b,o  +  b„q'  +  b,jq-^ 

Bjfq')  =  b^o  +  bi,q '  +  bj2q-2 
y(k)  is  the  excess  oxygen  level 

u,(k)  and  UjCk)  are  die  fuel  flow  and  air  flow  commands,  respectively 
k  is  the  sample  number 
delay  n,=15  seconds 
n2=  10  seconds. 

The  initial  parameter  values  used  are: 

a,  =  -0.562 
32=0.207 
83= -0.0545 

b, o=  0.140 
b„  =  0.0233 
b,2= -0.00723 
b2o=  -0.1% 
bj,  =  -0.0659 
b22=  0.0183 

3.  Model  for  water  level  locq): 

A(q  ')y(k)  =  q-" 'B,(q'')u,(k) +q'"^B2(q  ')u2(k)  IBq  5] 


where 

A(q*')  =  1  +  a,q'  +  a2q‘*+  ajq'^ 

B,(q ')  =  b,o  +  b„q-‘  +  b,jq-*  +  b.^q^ 

Bjfq')  =  b2o  +  b2,q'  +  b22q-*+  bjjq'^ 
y(k)  is  the  drum  water  level 

u,(k)  and  U2(k)  are  the  observed  steam  flow  and  feed  water  flow  command 
k  is  the  sample  number  and  delay  n,  =  20  seconds 
nj  s  20  seconds. 
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The  initial  parameter  values  used  are: 


a,  *-1.3199 
a*  *0.5725 

a, * -0.2518 

b, o*  0.001793 
b„*  0.001406 
b,i* -0.001742 
b,3* -0.002207 
ba,= -0.00052 
bji*  0.001479 
b22= -0.000424 
b23*  0.000194 


3  PROTOTYPE  ADAPTIVE  CONTROL  SYSTEM 


Hardware  Description 

Microprocessor-based  control  systems  are  increasingly  used  in  industrial  processes.  A  wide  variety 
of  coiiq)uterized  systems  packaged  in  different  forms  are  available  for  central  heat  plant  controls. 
Selecting  a  ^stem  can  be  a  trying  process  because  of  the  many  vendors  available  and  the  rapid  pace  at 
which  digital  hardware  and  software  conqionents  are  changing.  In  addition,  determining  the  control 
features  and  performance  objectives  requires  a  thorough  analysis  based  on  technical  and  economical 
considerations.  After  much  research,  a  personal  con^uter  (PC)  based  system  was  chosen  based  on  the 
system’s  flexibility  and  cost  (Wohadlo  et  al.  1990).  This  system  uses  a  PC  with  thirdi)aity  control 
software  combined  with  front-end  input/output  (I/O)  signal  conditioning  to  perform  analog-to-digital  and 
digital-to-analog  conversions.  System  functions  include  data  collection,  report  generation,  PID/advanced 
control,  alarming,  and  operator  interface.  The  control  software  selected,  THE  FIX^  is  produced  by 
Intellution,  Inc.  The  controller  used,  (MAC,  was  manufactured  by  Analog  Devices,  Inc.* 

THE  FIX  software  enables  the  fiMAC  controller  to  monitor  and  control  industrial  processes  when 
used  with  a  DOS-compatible  PC.  THE  FIX  is  a  menu-driven,  multitasking  software  package  that 
communicates  widi  distributed  I/O  devices  and  provides  users  with  the  ability  to  use  a  PC  workstation  as 
a  siq>ervi8ory  host  computer.  The  base  configuration  includes  software  for  analog  and  digital  I/O,  alarm 
detection  and  messaging,  real-time  trending,  online  calculation  blocks  and  user-drawn  color  graphics.  THE 
FIX  is  also  available  in  runtime-only  systems  that  are  useful  for  multiple  installations. 

The  pMAC-6000  is  a  modular  I/O  processor  that  offers  input  and  output  interfeces,  and 
communications  and  conqniting  ciqwbilities.  This  processor  can  be  us^  as  a  standalone  system,  with  a 
host,  or  as  a  part  of  a  distributed  control  system.  The  base  system  has  24  analog  and  48  digital  I/O  points, 
and  can  be  expanded  to  over  200  analog  and  1000  digital  I/O  points  per  cluster.  The  controller  central 
processing  unit  (CPU)  is  based  on  Intel  80188  microprocessor  ruiming  at  a  clock  rate  of  8  MHz.  The 
CPU  comes  with  2S6K  bytes  of  CMOS  battery-backed  RAM  memory  and  192K  bytes  of  ROM  memory. 
In  the  C  Programmable  version,  64K  bytes  are  taken  up  by  the  compiler  and  128K  is  left  for  the  user’s 
compiled  C  program.  Application  development  on  the  pMAC-6000  is  simplified  by  the  use  of  high-level 
languages,  special  libraries  of  high-level  functions,  program  development  aides,  and  runtime  communica¬ 
tion  stqqmrt  MCConun->-  is  die  Analog  Device’s  communication  protocol,  created  to  provide  runtime 
su^xirt  from  host  computers  to  pMAC  systems.  Intemqits  provide  a  mechanism  for  the  controller  to  react 
to  events  haiqiening  in  real  time  and  to  help  provide  structure  in  programs  designed  for  control  of  events. 
Six  interrupts  are  built  in  the  pMAC,  all  of  which  are  available  for  program  use. 

For  a  prototype  control  system,  issues  related  to  long-term  safety  and  reliability  needed  to  be 
addressed.  Loss  of  power  and  software  upset  could  cause  the  command  signals  (water  flow,  air  flow,  and 
fuel  flow)  to  be  turned  off,  either  singly  or  as  a  group.  In  particular,  the  software  accomplishes  fuel/air 
intalocking.  A  software  crash  could  cause  a  loss  of  air  flow  while  fuel  flow  continued.  A  backup  system 
to  address  these  problems  had  to  be  inqilemented  before  testing.  Redundant  power  supplies  were  needed 
to  provide  DC  power  to  the  backiq)  systeriL  A  7312  LOOPMATE  manufactured  by  Control  Technology 
Inc.  (5734  Middl^rook  Pike,  Knoxville,  TN,  tel  615/584-0440)  was  chosen  as  system  backup.  The 
LOOPMATE  is  designed  to  provide  automatic  and  hard  manual  backup  for  process  loop  control  in  harsh 
environments.  The  LOOPMATE  provides  a  highly  visible  display  of  critical  sensor  and  control  signal 


btfethaion,  Inc.,  315  Norwood  Park  South,  Norwood,  MA  02060,  tel  617/769-8878;  Analog  Devices,  Inc.,  One  Technology 
Way,  PO  Box  9106,  Norwood,  MA  02060,  tel.  617/326-6666.  Note  diat  uMac  controllers  are  currently  sold  by  Azonix  Corp., 
900  Middlesex  Turnpike,  Building  G,  Billercia,  MA  01821,  tel.  1-800«36S-1663  (toll-free). 
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variables  and  allows  (^onal  remote  set-point  changes  and  display.  During  nonnal  operation,  the 
LOOPMATE  is  transparent  to  the  controller.  The  LOOPMATE  takes  over  wtren  an  out-of-control 
situation  is  detected.  The  user  may  then  select  one  of  the  three  manual  control  modes,  either  to  default 
to  the  last  valid  control  variable  ouqiHit,  ramp  to  a  i»eset  control  variable  output  value,  or  default  to  4  mA 
or  1  V. 

To  test  the  prototype  control  system  in  an  opertding  boiler,  it  was  first  necessary  to  obtain  detailed 
existing  plant  pneumatic  control  system  data,  including  construction  data  and  schematics,  ladder  diagrams 
of  electrical  safeties,  plant  control  requirements,  and  control  system  characteristics.  The  pMAC  controller 
needs  to  have  current  inputs,  so  pneumatic-to-current  converters  had  to  be  installed  in  the  control  lines. 
The  oxygen  level  signal  was  already  a  1  to  S  V  signal,  but  the  other  signals  were  3  to  IS  psi  pneumatic 
signals.  Further  processing  was  necessary  on  air  flow,  feedwater  flow,  and  steam  flow  signals.  These 
signals  were  measured  with  differential  pressure  flowmeters,  so  a  root  extractor  was  necessary  to  obtain 
the  true  flow.  Control  was  switched  over  from  the  existing  control  system  to  the  pMAC  by  a  multipon 
valve  installed  by  a  plant  technician.  To  switch  control  over  horn  the  automatic  coittrol  to  the  pMAC, 
the  control  system  was  put  in  manual  mode  and  the  valve  was  switched.  When  the  pMAC  signal  and  the 
manual  command  signal  matched  to  provide  a  bumpless  transfer,  the  panel  was  switched  back  into  the 
automatic  mode,  allowing  the  pMAC  to  take  over  system  control. 


GPC  Algorithm 

Basic  Controller  Structure 

The  purpose  of  any  controller  is  to  ensure  robust  performance.  This  means  that  the  controller  must 
provide:  (1)  good  tracking  of  the  reference  input  by  f^ant  output,  and  (2)  disturbance  rejection,  i.e.,  to 
maintain  tracked  values  in  spite  of  external  dishirbances  and  to  preserve  such  performances  under 
variations  in  plant  characteristics.  For  a  large  class  of  processes,  this  purpose  has  been  demonstrated  to 
be  well  served  by  robust  self-tuning  controllers.  A  self-tuning  controller  (Figure  1)  consists  of  the 
following  blocks: 


Figure  1.  Diagram  of  a  Self-Tuning  Controller. 
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1 .  Regulator/tracker.  This  block  is  an  algorithm  that  computes  corttrol  signal  on  the  basis  of  (a) 
feedback  signal  (feedback  sensor  output),  (b)  reference  input  (preplanned  set  point  or  time  varying  signal, 
or  a  priori  unknown  input  signal),  and  (c)  feedforward  signal  (feedforward  sensor  output). 

2.  Control  design  block.  This  block  is  conceptual.  It  includes  (a)  on-line  regulator/tracker 
reparameterization  that  is  taking  place  at  every  sampling  period  on  die  basis  of  current  plant  parameters 
cmnputed  by  parameter  estimator,  and;  (b)  an  off-line  regulator/tracker  design  using  an  underlying 
controller  synlfaesis  method  such  as  GPC  and  setting  controller  features,  such  as  regulator  structure, 
prediction  atul  control  horizons,  and  rate  limits  on  control  signal.  However,  if  hardware  permits, 
regulator/tracker  redesign  can  be  done  in  real  time  in  closed  loop  to  allow  the  on-line  change  of  controller 
feature. 

3.  Parameter  estimator.  Within  a  specified  model  structure,  this  algorithm  computes  the  parameters 
of  the  model  to  approximate  plant  input/output  behavior  on  the  basis  of  output  measurements  (feedback 
ouQnit  sensor  sig^)  and  plant  input  (control  signal). 

4.  Siqrervisor.  This  block  includes:  (a)  the  real-time  diagnostics  of  the  self-tuner,  such  as  identifier 
conditioning,  depoiding  on  the  richness  of  online  information  in  the  input/output  signals,  (b)  non-real-time 
operator/designer  interfoce  with  an  identifier  for  changing  identification  procedure  and/or  identifier 
structure,  (c)  non-real-time  operator/designer  interfoce  witii  control  design  block  for  resetting  ^iuning 
knobs”  in  control  design  procedure  and  the  regulator/tracker  structure.  If  hardware  permits,  tiie  last  two 
tasks  can  be  performed  in  real-time,  resulting  in  an  “intelligent”  self-tuner. 


The  GPC  Methodology 
Design  Components 

The  synthesis  of  the  GPC  regulator  involves  the  following  components: 

1.  A  model  that  expresses  the  iiqiut/ou^ut  plant  behavior  and  typical  plant  disturbances  in  a  finite 
set  of  parameters.  The  most  suitable  model  structure  for  a  large  variety  of  plants  is  of  the  form: 

y(t)+a,y(l-l)+-+a„y(t-n^  «b„u(t-l)+-.. 

+b^u(t-n^-l)  +a(t)  +c,  a(-l)  +-  +c„  a(t-n^) 

and 

a(i)-a(i-l)  =  ^(i)  [Eq  7] 


Where  u(i),  y(i),  and  a(i)  are  the  values  of  plant  input,  ou^ut  and  disturbance  at  instant  i,  respectively. 
Swpience  {^i)}“^  is  an  uncorrelated  zero  mean  random  pulse  sequence  with  variance  o^.  Consequentiy 
from  Equation  7  it  is  seen  tiiat  {a(i)}'°^  obtained  by  discrete  integration  from  {^i)}*H»  is  a  stq) 
sequence  suitable  for  nuxleling  the  step-like  disturbances.  By  definiqg  q  as  a  forward  time  shift  operator 

qy(t)=y(t  +  l)  [Eq  8] 
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aiMl  q  '  as  a  backward  time-shift  operator: 


q-ly{t)  =  y(t-l) 

mid  then  by  combining  Equations  6  and  7  and  factoring  out  y(t)  and  u(t-l),  the  so  called  Controlled  Auto- 
Regressive  Integrated  Moving  Average  (CARIMA)  model  is  obtained; 

A(q-')y(t)  =  B(q-')u(t-l)  ♦CCq  ')!^  [Eq  10] 

A 


where: 


A{g-')  «  1 +fl,q 

A  =  l-q~\difference  operator). 

The  model  parameters  are  updated  online  by  a  parameter  estimator  that  uses  an  a  priori  fixed  model 
structure.  The  model  structure  (orders  of  polynomials  A(q''),  B(q''),  and  C(q''))  can  be  changed  on-line 
by  the  siqiervisor  block  or  offline  by  the  designer.  The  model  primarily  serves  as  a  basis  for  the  optimal 
prediction  at  time  t  of  the  future  plant  outputs  y  (t  j),  1,2, ...,  to  the  present  and  future  controls. 

2.  A  future  reference  sequence.  {W  (t  +  j)},  J  -  I,  2, ....  This  might  be  either  a  fixed  setpoint  or 
ptcphmned  sequmce.  However,  the  design  procedure  can  be  augmented  to  design  the  regulator/tracker 
for  an  a  priori  unknown  reference  sequence,  as  well. 

3.  A  measure  of  the  cost  of  the  future  predicted  output  errors  {e  (t  +  j)}  and  the  control  effort 
{u  (t  +  j  -  1)},  is  usually  provided  by  a  functional  J  ({6  (t  +  j)},  {u  (t  +  j  -  I)}).  This  measure  denotes 
the  optimal  predicted  value  of  y  (t  +  j)  given  data  iq>  to  time  t  by  y  (t  j  1 1),  and  the  future  control 
incremoits  (control  “moves”)  by  au  (t  +  j)  =  u  (t  +  j)  -  u  (t  +  j  -  1).  In  the  GPC  methodology,  such  a 
function,  often  sinqily  referred  to  as  a  cost-function,  is  given  by: 

N 

JoPc(NpN2,N„,>.)-f;  [y(t+j|t)-W(t+j)rU^[Au(t+j-l)]^  11] 

j-N,  j»1 

Where 

e  (t  +  j)  =  y  (t  +  j  1 1)  -  W  (t  +  j)  is  the  predicted  output  error  at  instant  t+j 

N,  is  the  cost-starting  error  horizon 

Nj  is  tile  ouquit  prediction  horizon 

Ny  is  the  costing  control  horizon 

X,  is  the  control  weighting. 

The  assignment  of  costing  hwizons  (the  boundaries  of  the  time  intervals  where  the  cost  is  computed) 
for  the  output  error  and  the  input  signal  according  to  Equation  1 1  is  done  either  offline  1^  the  designer, 
or  online  a  supervisor  block.  The  quadruple  [N,,  N2,  Ny,  A.]  can  be  viewed  as  a  set  of  “tuning  knobs” 
for  a  GPC  algorhhm. 
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4.  An  optimization  procedure  that  at  every  current  time  moment  t  selects  the  future  control 
sequence  {u  (t  j  -  1)},  j  =  I,  2,  ....  to  minimize  the  cost  function  Jqpc-  Control  sequence  can  be 
determined,  therefore,  not  further  than  the  output  prediction  horizon,  i.e.,  max  {j}  =  Nj.  However,  a 
simplified  assumption  is  often  made  after  the  (t  N„  —  1)  steps-control  increments  are  zero,  which  is 
equivalent  to  assigning  an  infinite  cost  to  control  “moves”  after  the  (t  N„  -  1 )  stq>s.  If  an  optimization 
procedure  is  unconstrained,  then  it  is  possible  to  derive  an  explicit  expression  for  the  analytical  calculation 
of  the  optimal  control  sequence,  lliis  allows  a  real-time  implementation  of  GPC.  Explicit  expression  can 
also  be  drived  for  some  cases  of  constrained  optimization,  such  as  the  rate  limitation  on  the  control 
signal,  using  the  Lagrange  multiplier  techniques  (Wilde  and  Beighter  1967). 

5 .  Signal  conditioning  filters  increase  robustness  of  the  algorithm  against  unmodeled  dynamics  and 
disturbances.  In  ftict,  the  polynomial  C(q'')  in  the  plant  model  can  be  viewed  not  as  a  disturbance  model, 
but  as  a  filter  (often  noted  as  T(q'')  polynomial)  that  must  have  roll-off  at  high  frequencies  to  reduce 
predicted  error  due  to  unmodeled  dynamics  that  mi^t  not  be  captured  well  at  high  frequencies  by 
Controlled  Auto-R^ressive  and  Integrated  Moving  Average  (CARIMA)  models.  Also,  A  can  be  viewed 
as  an  integrator  that  eliminates  prediction  error  due  to  the  model’s  inaccurate  DC  gain,  and  removes  the 
effect  of  DC  load  disturbance  on  the  output  error. 

Formulation  of  the  Optimal  Predictor 

The  key  element  in  the  GPC  strategy  is  the  optimal  j-step-ahead  predictor,  an  algorithm  that  uses 
all  available  input/output  measurements  up  to  and  including  time  t,  {y  (t),  y  (t  -  1), ...,  y  (0),  A  u  (t  -  1), 
A  u  (t  -  2), ...,  A  u  (0)},  as  well  as  the  future  control  increments,  {Au  (t),  Au  (t  +  1), ...,  A  u  (t  +  j  -  1)}, 
for  computing  die  future  output  estimate  y  (t  '**  j  1 1),  which  minimizes  the  variance  of  the  j-step-ahcad 
prediction  error,  e,(t  +  j  1 1)  *  y  ( t  +  j)  -  S' (t  +  j  1 1),  where  y  (t  +  j)  is  the  true  future  j-step-ahead 
output 

Below,  the  j-stqi-ahead  predictor  is  derived  for  the  simplified  case  of  nc=0,  i.e.,  C(q'')=l.  In  this 
case,  the  CARIMA  model  takes  the  form 

A(q-')y(t)  =B(q->)u(t-l)  *1^(1)  [Eq  12] 

A 

The  derivation  for  n(>0  is  given  in  Favier  (1987).  Following  Ljung  (1991),  the  Dit^handne  identity: 

1  =E.(q-')A(q-')A  ^q^F;(q-')  [Eq  13] 


is  solved  for  Ej(q')  and  Fj'(q''),  where  Ej(q  ')  and  Fj'(q  ')  are  polynomials  of  degree  j-1  and  n,-l, 
respectively,  which  are  unique^  defined  given  A(q'')  and  the  prediction  interval  j.  Multiplying  Equation 
12  by  q*  Ej  (q  ')  A  yields 

Ej(q’')A(q-')Ay(t+j)  »  E.(q-')B(q'')Au(t+j-l)  +E.(q-')4(t+j)  [Eq  14] 

Rewriting  Equation  13  as 


E3(q-•)A(q-')A  =  I-q^F;(q-') 


[Eql5] 
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and  substituting  Equation  IS  into  the  left  hand  side  of  Equation  14  yields: 


[l-q^F;(q-')]  =  Ej(q-‘)B(q-')Au(t+j-l)+E.(q-‘)^(t*j)  [Eq  16] 

Defining  Gj(q'')sEi(q'')B(q''),  Equation  16  can  be  rewritten  as 

y(t*j)  =G.(q-')Au(t*j-I)  ♦F;(q-')y(t)  ^E.(q-')4(t+j)  [Ed  1^] 

Substituting  q=l  into  Equation  13  yields  Fj'=l,  Fj'(q’'),  which  can  be  factored  as 

F/Cq-')  =  l+F/q-')A  [Eq  18] 

where  F|'(q'')  is  a  potynomial  in  q''.  Sidistituting  Ec^tion  18  into  Equation  17  yields 

y(t+j)  »  G.(q'')Au(t+j-l)  +y(t)  +Fj(q-’)Ay(t)  +E(q-')^(t+j)  [Eq  19] 


Since  Ej  (q  ')  is  a  polynomial  of  degree  j-l,  the  last  term  in  Equation  19  is  unpredictable;  therefore  the  best 
estimate  y  (t  +  j  1 1)  of  y  (t  +  j),  given  current  and  past  output  and  control  values  as  well  as  future  control 
increments  Au  (t  j  —  1),  i  >  1,  is  obtained  by  dropping  fte  last  term  in  Equation  19: 

y(t+j  |t)  =  G.(q-')Au(t+j-l)  +y(t)  +F.(q-')Ay(t)  [Eq  20] 

This  is  the  so-called  incremental  predictor,  which  predicts  the  fiill  value  of  y(t^j)  using  the  input  and 
output  increments,  Au(i)  and  Ay(i).  Denoting  the  coefficients  of  polynomials  Gj(q'')  and  Fj(q'')  as  gj^, 
i=0,l, ...,  iiy  and  i=0,  1, ...,  n,-l,  respectively.  Equation  20  can  be  rewritten  as: 

y(t+j  |t)  =  Au(t+j-i)+p^  [Eq  21] 


[Eq23] 


which  corresponds  to  the  predicted  free  future  response;  and  (2)  the  predicted  forced  future  response: 

j 

which  is  the  convolution  of  the  future  control  incronrats  with  plant  step  response,  since  the  first  j 
coefficients  of  G/q  *)  are  just  ordinates  of  the  plain  step  response.  Combining  all  the  predictions 
{^(t+j  1 1),  j=  1,  ....  N,}  into  a  vector  f  obtains  the  key  equation: 

Y  -  Gu  +  P 


where 

t  -  mt+i|t),  ^(t+2|t),  •••,  j(t+Njt)r 
0  «  [Au(t),  Au(t+1),  Au(t  +  Nj-l)r 
P  •  tPp  Pj.  •  Pn,]^ 

■g.  0  ...  O' 

^  82  gi  ■"  0 

u  * 

•M  «««  «««  «•« 

8n,  Sn,  "■  81 
L  '  ’  J 


Equation  24  is  the  basis  for  the  dwivation  of  the  GPC  control  algorithm. 


Formulation  of  the  GPC  Control  Law 

Figure  2  shows  the  conceptual  framework  for  the  derivation  of  the  GPC  control  law.  The 
derivation  starts  with  obtaining  a  sequence  of  future  control  incranents  Au(t+j-l),  j=  1,  ...,  N^,  which 
minimize  a  cost  function  (Equation  1 1).  By  taking  into  account  costing  GPC  horizons,  the  compon«its 
of  Equation  24  are  given  by: 


t  -  mt+NJt),  Kt+N,  +  llt),  -,  ^(t+NjOr 


[Au(t), 

Au(t 

^1),  •, 

Au(t  + 

[Pn,»  Pn 

-i*i» 

PN.r 

8n, 

Q 

»N,-2 

0 

8n+1 

«N. 

8n.-, 

0 

8n,-i 

8n,-3 

8n,-3  ■“ 

8n,-n 

8n,-i 

8n.-2  ••• 

8n,-n. 
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Figure  2.  A  Conceptual  Framework  of  the  GPC  Law. 


Introducing  vector 

W  =  [W(t+N,),  W(t+N,  +1),  •••,  W(t-*-Nj)r  [Eq  25] 


Equatioi  1 1  can  be  rewritten  in  vector-matrix  form  as 


Jqpc  =  (W-t)T(W-t)  +>.0^0=  (W-GU-Pf  (W-GU-P)  A,U  ^0  l^q  26] 


Ugpt,  which  provides  the  global  minimum  of  quadratic  function  Equation  26,  is  found  by  solving: 

=0  [Eq  27] 

ao 


which  yields  a  computation  formula  for  a  single  step  of  the  GPC  control  law: 

0^  =  (G  ^G  +  XT)-'  G  T(W-P)  nSq  28] 

The  GPC  control  law  is  generated  using  the  “receding  horizon  strategy.”  At  each  sample-instant  t, 

1 .  Updated  model  parameters  are  passed  from  parameter  estimator  to  update  the  parameters  of  the 
(nedictor  and  the  single-step  control  algorithm.  (In  the  nm-self-tuning  GPC  version,  this  step  is  omitted.) 

2.  From  the  current  measured  output  y(t)  and  die  previous  values  of  inputs  and  ouqmts,  vector  P 
of  die  predicted  free  plant  response  is  computed  and  vector  W  of  the  iiiture  reference  input  sequoice  is 
evaluated 


3.  The  future  optimal  control  vector  is  computed 

4.  Only  the  first  element  Au(t)  of  vector  U,p,  is  used  to  generate  and  ^ply  the  actual  control  signal 
u(t)=u(t*l)+A(t) 

5.  At  the  next  time  moment,  ta,j+l=t^,  the  whole  procedure  is  repeated  to  account  for  a  new 
output  measurement  y(t^)  and  to  compute  u(t^)=u(tnew-l)+Au(t^). 

Presently,  efficient  recursive  algorithms  can  compute  P  via  the  multistep  predictor  formulation  and 
the  inversion  of  matrix  Q*G^G+XI  (Favier  1987).  Multivariable  GPC  algorithms  have  also  been 
developed  (Lambert  1987;  Favier  1987). 

Recursive  Least  Squares  (RLS) 

Considering  CARIMA  model  (Equation  10),  which  can  be  rewritten  as  an  incremental  model: 

A(q-')Ay(t)  =  B(q-')Au(t)  ■^C(q-')^(t)  t^q  29] 


where: 

Ay(t)  =  y(t)-y(t-l) 
Au(t-k)  =  u(t-k)-u(t-k-l) 


If  Equation  29  is  ^ut  into  a  form  that  will  fit  into  a  Recursive  Least  Square  (RLS)  (Favier  1987; 
Ljung  1987;  Ljung  and  Astrdm  1987): 


where: 


y(t)  =e'*(t)*4(t) 


[Eq30] 


0  (aj,  32,  ■*',  bj,  b|,  Cq,  C],  "•] 

composes  of  polynomials  A,  B,  and  C  coefficients; 

♦(0^-A></-l),-A></-2),-^Att(/-l),A«(/-2),...,M</-l),w(/-2),--] 
con^ses  of  input-output  measurements; 
w(t)  is  an  estimate  of 

RLS  is  based  on  the  predicted  error  method.  The  RLS  algorithm  is  the  result  of  the  minimization 
of  the  following  cost  function  with  respect  to  6  (Warwick  1988): 

J  =  l^X,'-'e*(i)  [Eq31] 
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where 

e(t)  is  the  estimated  error  between  the  current  output 
(y(t))  and  the  estimated  output  (^(t)) 

X  is  the  forgetting  factor. 

Equation  31  can  be  rearranged  as: 


j  =  4rx-  *(y(t)-(?^^(t))^ 


[Eq  321 


The  forgoing  factor  X  is  a  data  weighting  factor  that  weights  the  current  data  more  than  the  old 
data.  That  is  to  say,  it  forgets  the  old  data.  This  is  an  ad-hoc  method  for  dealing  with  uncertain 
estimation.  Finally,  minimization  of  Equation  33  gives  (Warwick  1988): 


d(t)  =  d(t-l)+P(t)0(t)(Ay(t)-d'(t-l)0(t)) 
p-‘(t)  =  p-‘(t-l)+^(t)^'(t) 


[Eq  33] 
[Eq34] 


where,  P(t)  is  the  covariance  matrix.  The  covariance  matrix  P(t)  is  a  measure  of  the  convergence  of  the 
estimated  parameters.  In  general,  if  the  trace  of  P(t)  is  small,  generally  less  than  one,  the  estimated 
parameters  will  approximate  the  true  parameters.  If  the  trace  of  P(t)  is  large,  the  estimated  parametm 
are  uncertain. 

The  conq)utation  of  Equation  34  and  Equation  35  is  expensive,  since  at  each  sampling  interval,  a 
matrix  inversion  is  performed.  However,  the  matrix  inversion  lemma  presented  in  AstrOm  and 
Wittenmark+  (1989)  can  be  employed  to  reduce  the  computing  time  and  improve  nummcal  accuracy, 
from  the  lemma: 


(A  +BCD)*  =  A  ‘ -A  *8(0  *  +DA  >B)  *DA  > 


[Eq351 


where: 

A,  C,  and  (C‘  +D  A  '  B)  are  nonsingular  square  matrices.  That  is,  inversion  is  substituted  by 
multiplication  and  addition.  The  matrix  inversion  lemma  is  next  applied  to  Equation  35.  This  leads  to 
the  standard  RLS  (Warwick  1988): 

estimated  error, 


e(t)  =  Ay(t)-$(t-l)<A(t) 


Kalman  gain  update. 


K(t)  = . _ 

[X+^'(t)P(t-l)«(t)I 


[Eq36] 


(Eq371 
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Covariance  matrix  iqxlate. 


P(t)  »  0^(t)»'(t)P(t  l)j  1  Pq  3gj 

X>*'(t)P(t-l)*(t)  X 


estimated  parameters  update. 


6(t)  -  6(t-l)  +K(t)e(t) 


[Eq39] 


That  is. 


[new  estimates]  =  [old  estimates]  *  Kalman  gain  *  [estimated  errors] 


[Eq40] 


Equation  38  is  not  numerically  well-conditioned.  A  better  way  to  update  the  covariance  matrix  P(t) 
is  to  iqrdate  the  square  root  of  P  instead  of  P.  Rewrite  P  >  UDU’,  where  U  is  the  upper  triangular  matrix 
of  P  and  D  is  the  diagonal  matrix  of  P.  The  square  root  of  P  becomes  UD'*^.  This  procedure  is  referred 
to  as  UD  factorization  (Bierman  1977). 

In  general,  the  RLS  algorithm  does  not  estimate  the  C  (color  noise)  polynomial.  An  extension  of 
the  RLS  algorithm  will  approximate  the  random  noise  term.  This  is  called  extended  RLS  (ELS)  (Ljung 
and  Astrdm  1987).  It  is  irtqiossible  to  estimate  iwdom  white  or  color  noise  directly.  However, 
introducing  another  variable,  w,  obtains  (Ljung  and  AstrOm  1987): 

w(t)  =  Ay(t)-d'(t-I)*  (t)  ^ 

Equation  38  will  give  the  approximate  noise  at  t.  That  is,  w(t)  will  be  approximately  equal  to  ^(t).  For 
sinq>licity,  RLS  and  ELS  will  be  considered  synonymous  for  the  balance  of  this  report. 

The  RLS  algorithm  caimot  track  slow  time-varying  parameters  because  when  the  estimated  error 
becomes  small,  so  do  the  covariance  matrix  P  and  Kalman  gain  K.  This  is  one  of  the  major  drawbacks 
of  RLS,  which  does  not  retain  alertness  or  adaptivity.  When  the  parameters  drift  (i.e.,  when  the  estimated 
error  is  large),  the  P  matrix  and  K  do  not  change.  Therefore,  the  RLS  algorithm  gives  poor  and 
sometimes  wrong  estimates  for  time-varying  systems.  However,  there  is  a  parameter  available  in  the  RLS 
algorithm  that  can  be  used  to  address  tUs  problem — the  forgetting  fector  X.  Even  though  the  forgetting 
ftictor  can  take  on  any  value  between  zero  and  one,  it  has  been  found  that  0.90  can  be  considered  as  a 
lower  bound  and  that  0.95  is  a  good  starting  or  default  value  in  industrial  practice  (Warwick  1988).  In 
general,  the  constant  forgetting  ftictor  will  improve  the  estimation  of  slow  time-varying  systems.  However, 
at  the  steady  state  (i.e.,  where  the  parameters  do  not  vary  or  vary  very  slowly)  wi&  disturbances,  the  RLS 
algorithm  may  cause  the  estimated  parameters  to  oscillate.  Therefore,  the  RLS  with  constant  forgetting 
fimtor  is  not  recommended  for  long-term  identification. 

There  are  many  ad-hoc  procedures  to  keep  the  RLS  algorithm  “awake,”  including  leakage,  constant 
trace,  covariance  resetting,  dead  zone,  and  others.  Howev^,  these  ad-hoc  procedures  require  process 
knowledge  to  work.  For  example,  tire  RLS  with  covariance  resetting  prevents  the  covariance  P  from 
ftdling  below  a  certain  level.  If  the  covariance  P  is  below  that  level,  it  is  reset  to  some  large  value.  The 
lower  bound  of  P  and  the  value  to  which  P  is  reset  are  process-dependent  Human  operators  must  enter 
die  process  information  before  the  RLS  algorithm  starts.  Such  ad-hoc  procedures  are  not  very  effective 
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in  devel<^ing  a  general  puipose  identification  scheme.  Another  way  to  modify  the  RLS  algorithm  is  to 
allow  the  foigetting  factor  to  vary.  The  variable  forgetting  foctor  is  changed  according  to  the  estimated 
error.  If  the  estimated  error  is  large,  X(t)  will  be  set  to  a  small  value.  If  the  estimated  error  is  small,  A.(t) 
will  be  set  to  one.  By  rqrlacing  Equations  39  and  40,  the  RLS  with  variable  forgetting  fiu:tor  becomes 
(Warwick  1988): 


[y(t)-V(t)6(t-l)y  [Eq  42] 

[l^^'(t)P(tl)*(t)]o 

K(t)  = _ _  [Eq  43] 

[1 

W(t)  =  P(t-l)-K(t)#'(t)P(t-l)  [Eq44] 

where: 

X(t)  >=0.90  for  practical  implementation 

o/a««1000,  o„  is  the  variance  of  the  zero  mean  noise  sequence. 

The  above  development  assumed  an  input  signal  that  is  persistently  excited,  i.e.,  diat  the  irq}ut  signal 
contains  all  frequency  modes.  Under  feedback  control,  there  is  no  guarantee  that  the  ii^ut  signal  is 
persistently  excited.  When  the  input  signal  is  not  persistently  excited.  Equation  39  becomes  P(t)=P(t-l)/X. 
Therefore,  the  covariance  matrix  P  will  “blow  if  the  input  signal  is  not  “rich  enough.”  This 
phenomenon  is  usually  referred  to  as  estimator  “windup”  (Harris  and  Billings  1 98S).  A  detailed  discussion 
will  be  presented  in  the  self-tuning  control  section  after  the  introduction  of  the  generalized  predictive 
control  (CPC)  law. 

Inclusion  of  Constraints  on  the  Control  Signal  in  the  GPC  Control  Law 

While  both  magnitude  and  rate  constraints  on  control  signal  might  be  needed  in  a  specific 
inq>lementation,  it  is  numerically  more  efficient  to  introduce  only  rate  constraints  that  usually  indirectly 
ensure  adequate  control  signal  magnitudes  as  well.  Thus,  it  is  assumed  diat  control  “moves”  must  stay 
within  the  prescribed  feasibility  bound  during  the  whole  costing  control  interval: 

ttj  ^  Au(t+j-l)  <  Pj,  j  =  1,  -,  N„  [Eq  45] 

To  derive  an  algorithm  for  the  computation  of  die  consumed  optimal  control  vector  denoted  as 
assume,  for  simplicity,  that  only  one  control  increment  of  U,  saturates,  i.e.,  it  needs  to  be  on  the  boundary, 
for  example: 

Au(t+j-l)=a.  [Eq46] 


In  this  case  it  is  possible  to  find  using  a  Lagrange  muhiplier  technique.  For  this  purpose, 
function: 


L  =  (G,U,+P-W)^(G,U,+P-W)  +A.0/U,  +2Mi[Au(t+j-l)-ap  Pq  47] 
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is  minimized  by  solving  d\Jd\i=^.  This  yields: 

C,  =  (G/G,  G,^(W-P)  -^(G  Jg,  -XI)-'  MjCj  [Eq  48] 

and 

Cj  -  [0,  1.  0.  0]^  [Eq  49] 

or 

Uc  =0^-(G,"G,-XI)-'M5e.  [Eq50] 

Combining  the  jth  equation  of  Equation  SO  with  Equation  46  yields 

tEq  51] 

where  h^  is  an  element  (j  j)  of  matrix  (G  jGi+XI)"'.  Equation  S2  is  solved  for  unknown  and  finally  C. 
is  obtained  from  Equation  SO.  This  procedure  generalizes  to  any  number  of  saturated  controls  in  vector 
tic  (Clarke  and  Tsang  1988).  A  heuristic  choice  of  saturated  elements  of  Ue  is  to  pick  those  elements  of 
that  violate  feasibility  boundaries. 

Commissioning  of  a  Self-Tuning  Regulator 

At  the  beginning  of  the  tuning  phase,  an  adaptive  control  law  should  not  be  used  because  the 
process  parameters  are  uncertain.  Therefore,  the  self-tuning  control  (STC)  requires  a  start-up  procedure. 

There  are  two  ways  to  start-up  the  STC:  open-loop  identification  and  closed-loop  identification 
(Mohtadi  1988).  For  the  open-loop  identification,  the  process  parameters  are  estimated  in  open-loop 
fashion  with  the  pseudo-random  binary  sequence  (PRBS)  as  the  input  signal.  In  general,  the  commission- 
h^  period  would  be  from  about  SO  to  100  samples.  However,  for  industrial  process  control,  die 
commissioning  period  may  require  several  hundred  samples  due  to  the  long  time  constants  of  die 
processes.  After  the  commissioning  period,  die  adaptive  controller  “kicks  in”  and  the  PRBS  is  cut  off. 

Closed-loop  identification  uses  a  proportional,  proportional-integral,  or  proportional-integral- 
derivative  controller  to  close  the  loop.  The  settings  for  the  controller  will  be  conservative.  However, 
pretuning  the  PID  controller  is  required  (Appendix  B).  Also,  an  additional  PRBS  signal  is  added  into  the 
control  signal  to  improve  the  estimation.  After  the  commissioning  period,  the  controller  and  the  PRBS 
are  disengaged.  If  the  process  is  open-loop  stable,  bodi  start-up  procedures  can  be  used.  However,  if  the 
process  is  open-loop  unstable,  the  closed-loop  identification  procedure  must  be  used.  After  the 
commissionirig  period,  die  input  signal  may  not  be  persistentiy  excited.  A  good  control  law  will  try  to 
make  the  control  signal  as  smooth  as  possible.  Unfortunately,  identification  needs  a  persistently  excited 
plant  to  perform  well.  At  the  steady  state,  a  probing  signal  (e.g.,  PRBS,  periodic  signal,  etc.)  is  added  to 
die  control  signal  to  ensure  the  performance  of  the  estimator  (McDonald,  Kwatney,  and  Spare  1971). 
However,  the  probing  signal  may  cause  the  output  to  ring.  If  the  oscillation  of  output  cannot  be  tolerated, 
the  estimator  must  be  turned  off.  STC  always  has  conflicting  objectives.  The  estimator  performs  well 
with  persistent  excitation,  but  good  control  dictate  keeping  the  input  signal  as  constant  as  possible. 
Therefore,  there  is  a  constant  trade-off  between  identification  and  control  requirements  and  objectives. 
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Program  Deflcriptknis 

The  Software  Structure  of  the  Control  Algorithm 

The  control  algorithm  is  a  timer-initiated  interrupt  procedure  operating  on  the  pMac  when  the  pMac 
is  in  control  mode.  The  pMac  interrapt  system  remains  enabled  during  most  of  the  algorithm  to  permit 
plant  data  scanning  to  proceed  (continued  scanning  will  permit  the  implementation  of  digital  filtering  of 
the  input  data  stream  at  some  future  time).  The  interrupt  system  is  dis^led  (»ily  while  current  plant  data 
are  cq;>ied  to  buffers  used  for  the  control  computation.  Copying  the  scanned  data  assures  that  die  plant 
data  used  in  the  control  computation  will  not  be  nKxiified  as  die  computation  inoceeds. 

The  sequence  of  the  control  algorithm  follows.  The  data  needed  for  the  steam  pressure  lo(^  is 
copted  first,  and  then  die  firing  rate  to  steam  pressure  controller  (GPC)  is  executed.  This  yields  a  firing 
rate  command,  used  to  calculate  fuel  flow  ba^  on  fuel  energy  contort  calculations.  The  resulting  fiiel 
flow  is  interlocked  against  available  airflow  (to  prevent  incomplete  combustion)  and  then  the  desired  fuel 
flow  command  is  sent  to  the  fuel  pump.  Next  current  plant  data  are  again  copied  to  the  control  buffer, 
and  a  ccuitrol  calculation  for  the  excess  oxygen  loop  is  performed.  This  calculation  remits  in  a  new 
command  signal  for  the  blower  dampers.  Finally,  current  plant  data  are  copied  one  final  time,  die  water 
level  control  is  computed,  and  the  water  pump  setting  comnuuid  is  issued. 

Figures  3  to  S  show  the  main  control  loops  (outer  loops)  fcH'  drum  pressure,  oxygNi  level,  and  drum 
water  level.  Each  of  these  main  loops  contains  additional  controllers  around  the  actuators  (fuel  flow,  air 
flow,  and  feedwater  flow  inner  control  loops)  to  assure  that  the  command  signal  sent  by  the  controller  is 
the  signal  sent  by  the  actuate.  These  inner  loops  are  displayed  in  Figures  6  to  8.  It  is  noted  that  all  duee 
main  control  loops  (drum  pressure,  excess  oxygen  level,  and  drum  water  level)  are  computed  using  the 
same  GPC  function  (widi  different  parameters,  of  course).  This  structure  permits  easy  substitution  of 
alternate  ccmtrol  strategies. 


Figure  3.  Drum  Pressure  Contrtd  Loop. 
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Figure  4.  Oxygen  Level  Control  Loop  With  Interlocks  and  Feedforward. 


Figure  5.  Water  Level  Control  Loop  With  Feedforward  Blocks. 


Controller  Structure  and  Logic 

Feed  Foiwaid  Loops.  For  a  perfectly  known  model  and  perfectly  measured  distuibance,  the  ou^t 
of  foe  fduit  can  be  made  precisely  equal  to  foe  setpoint  regardless  of  foe  type  of  distuibance,  wifoout  foe 
need  for  feedback  control.  However,  even  in  our  idealized  model,  there  are  odier  nonmeasuied 
distuibances.  For  examine,  coupling  between  loops  wUl  produce  an  effect  in  subsystems  indirectly 
affected  by  a  particular  ii^Mt  Specifically,  in  foe  nonlinear  model,  the  header  pressure  is  indirectly 
affected  by  changes  in  water  flow.  Anofoer  example  of  nonmeasuied  disturbance  is  noise  in  sensors  and 
actuators.  Becuise  of  foese  ncnuneasured  distuibiuices,  feedback  control  is  still  necessary.  The  use  of 
feedforward  imi»oves  the  regulation  since  there  is  no  wait  for  foe  error  to  show  up  in  the  output  before 
hddng  corrective  actions.  Also  foe  process  of  idontifying  plant  parameters  for  self  tuning  is  made  simpler 
by  foe  partial  cancellation  of  disturbances. 


Figure  8.  Feedwater  Flow  Inner  Control  Lo(h>> 


A  feedforwaid  controller  was  chosen  for  the  fuel  flow/header  pressure  loop  to  reject  the  effect  of 
a  varying  “steam  demand”  on  the  system.  Consider  steam  demand  as  the  position  of  an  imaginary  valve 
in  the  output  steam  line.  Steam  demand  is  one  of  the  inputs  to  the  nonlinear  model  of  Bell  and  Astrdm, 
and  is  retained  as  an  input  in  both  the  linear  coupled  and  decoupled  models.  The  steam  demand  is  not 
a  measurable  quantity,  and  dierefoie  cannot  be  us^  as  a  feedforward  disturbance.  The  outlet  steam  flow 
rale,  however,  can  be  measured,  and  the  transfer  functions  from  steam  demand  to  steam  flow  and  from 
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steam  demand  to  header  pressure  may  be  calculated  from  the  limar  model.  By  dividing  these  transfer 
functions,  the  overall  transfer  function  from  die  measurable  steam  flow  to  the  ou^t  header  pressure  is 
obtained.  This  transfer  fiuictian  is  used  in  the  feed  forward  controller  design  suggested  by  Stepham^lous 
(Wahlbeig  1990). 

Attempts  to  directly  t^ly  the  feedforward  design  techniques  described  above  to  the  water 
flowAvater  level  loop  could  possibly  fail  due  to  die  non-minimum  phase  (i.e.,  shrink  and  swell) 
characteristics  of  die  drum.  When  the  feed  forward  design  is  carried  out,  an  unstable  pole  tqipears  in  the 
feedforward  transfer  function.  Since  disturbance  signals  pass  throu^  this  transfer  fimction  in  the  c^n 
loop,  die  controller  cannot  be  strdiilized  by  die  plant  and  the  control  signal  saturates.  Alternative 
fee^orward  designs  for  controlling  drum  water  level  are  currendy  under  investigation. 

Tiiteriocks  and  Cross  Limit  Constraints.  To  «isure  that  excess  oxygen  is  always  maintained  at  or 
above  die  values  required  for  safe  operation,  and  if  possible  at  optimum  levels,  a  cross-limiting  inocedure 
was  implmiiented  to  calculate  fuel  flow  and  air  flow  control  inputs  to  the  plant  (Figure  9). 

The  excess  oxygen  seqwint  is  taken  as  the  optimum  excess  oxygen  for  the  current  fiiel  flow 
conditions  (Figure  10).  The  current  low  fuel  conditions  are  tditained  as  the  minimum  of  the  fuel  flow 
requested  for  tracking  the  header  pressure  and  the  fuel  flow  allowed  by  the  current  available  air.  The 
approximate  air  flow  needed  to  match  the  fuel  at  the  desired  firing  rate  is  calculated  from  combustion 
(fynamics  of  die  linear  model. 

The  constraint  on  fuel  flow  ensures  that  even  when  load  conditions  change  rapidly,  there  will  always 
be  sufficient  excess  oxygen  in  the  combustion  chamber  to  ensure  safe  operation  (if  necessary,  at  the 
expense  of  a  drc^  in  header  pressure).  Higher  efficiency  and  lower  pollutant  emission  are  achieved 
through  the  prc^  selection  of  the  excess  oxygtm  set  points,  which  are  fimctions  of  the  given  boiler  load. 

CXie  may  question  the  need  for  feedback  on  die  excess  oxygen  loq>.  It  may  appear  that  if  the 
interlocks  are  properly  designed,  the  current  air  flow  for  obtaining  the  optimal  excess  oxygen  under  the 
current  fuel  flow  conditions  will  be  asserted.  The  jmiblem  is  that  the  interlocks  rely  on  a  linear 
combustitm  model  that  may  not  be  valid  outside  a  certain  operating  range.  The  purpose  of  feedback,  then, 
is  to  eliminate  any  steady-state  error  in  excess  oxygen  due  to  nonlinearities.  However,  the  addition  of 


AF- Air  Flow  c  -  ControUn  value 

FF  -  Fuel  Flow  i  -  Interlock  value 

Figure  9.  Interlock  (Cross  Limit  Constraint)  Confi^iuration. 
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Figure  10.  Optimal  Excess  Oxygen  for  a  Given  Load. 

is  to  eliminate  any  steady-state  error  in  excess  oxygen  due  to  nonlinearities.  However,  die  addition  of 
feedback  on  the  excess  oxygen  loop  may  also  introduce  dangerous  transients  into  the  air  flow  signal, 
defeating  the  purpose  of  the  interlocks. 

Probably  the  best  solution  to  this  problem  is  to  use  feedback  with  the  interlocks,  but  to  place 
constraints  or  heavy  weights  on  the  control  signal.  This  may  be  performed  using  either  a  GPC  or  PID 
amtroller  and  should  eliminate  both  the  steady-state  error  and  the  risk  of  bad  transients.  Anodier  solution 
is  to  use  recursive  least  squares  identification  to  provide  the  interlocks  with  a  more  reliaUe  combustion 
model.  However,  this  approach  may  place  excessive  reliance  on  the  performance  of  the  identification. 
Further  testing  widi  various  boiler  models  will  be  necessary  before  a  definite  cotKiusion  can  be  reached. 

Control  System  Flow  Chart 

Figure  11  ^ws  die  main  program  flow  chart  Figure  12  shows  the  flow  chart  for  die  pMAC-6000 
ctmtroller  data-processing  procedure,  and  Figures  13  and  14  show  the  flow  charts  for  iruier  and  outer 
loops,  respectively. 


Startiqi  and  Shutdown  Procedures 

The  autmnadc  startup  and  shutdown  procedures  were  implemented  using  the  FIX  DMACS  software 
padcage  running  on  a  386SX  PC.  FIX  DMACS*  is  multitasking  software  made  up  of  foreground  tasks 
and  background  tasks  dud  provide  complete  monitoring  and  control  over  automated  processes.  The 
for^round  tasks  are  programs  that  operators  interact  widi  to  set  up  the  system,  display  d^  and  perform 
other  functions.  The  badtgtound  tasks,  such  as  Schedule  Reports  and  Tasks,  Alarm  Monitoring  and 
Processing,  once  started,  tun  without  intervention  from  the  operator.  The  FIX  DMACS  padcage  has  beoi 
installed  an  die  mieter  ctnittoi  console  of  the  prototype  system.  The  database  code  for  automatic  startup 
and  dwtdown  (Mocedutes  are  presented  in  Appendix  C. 


‘FIX  rafen  to  Fully  Inlegraied  Condpol  System:  DMACS  refers  to  Distributed  Manufacturing,  Automation,  and  Control  Software. 
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F^pre  12.  Flow  Chart  for  the  MAC  Data  Processiiig  Procedure. 
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Figure  13.  The  Flow  Chart  of  the  Inner  Lo(^  Control  Procedure. 


36 


Figure  14.  The  Flow  Chart  of  the  Outer  Loop  Control  Procedure. 
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Figure  14.  (Cont’d) 
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The  startup  and  shutdown  procedures  detailed  below  were  developed  in  accordance  with  NFPA  85B 
“Standard  for  Explosions  in  Natural  Gas-Fired  Multiple  Burner  Boiler-Furnaces”  (National  Fire  Protection 
Association  1989).  Observations  made  at  the  Abbon  Power  Plant  were  also  considered. 

Startup  Procedure 

The  startup  sequence  will  comply  with  NFPA  8SB  section  S-2.1.2,  and  will  have  to  be  customized 
for  every  plant  based  on  boiler  manufacturer’s  recommendations.  The  sequence  is  implemented  using  the 
FIX  DMACS  software,  and  can  easily  be  teconfiguied  to  suit  any  boiler.  A  generalized  startup  sequence 
is: 


1.  Conduct  a  prestartup  check  per  NFPA  8SB  section  5-2. 1.1.  In  particular,  bring  the  drum  water 
level  to  the  normal  operating  point  and  engage  the  drum  water  level  controller.  Open  dampers  and  verify 
an  open  flow  path  horn  air  irilets  to  the  stack.  Prompt  operator  if  manual  checks  are  required. 

2.  Start  fans  and  achieve  a  purge  air  flow  rate  in  accordance  with  NFPA  8SG.  Engage  and  use 
the  air  flow  control  loop  to  establish  the  air  flow  rate. 

3.  If  the  boiler  is  equipped  with  regenerative  type  air  heaters  and  gas  recirculating  fans,  start  them, 
or  prompt  the  operator  to  start  them. 

4.  Precharge  the  burner  header  with  gas  per  NFPA  8SB  section  5-2. 1 .2  item  d  and  check  for  leaks. 

5.  Perform  a  unit  purge. 

6.  Close  the  main  fuel  control  valve,  then  open  the  main  safety  shutoff  valve. 

7.  Establish  correct  header  pressure  for  burner  lightoff. 

8.  Establish  correct  igniter  header  pressure  for  ignition. 

9.  Adjust  die  air  register  or  damper  on  the  burner  selected  for  light-off  to  the  correct  light-off 
position. 

10.  AttemiA  to  ignite  the  first  igniter.  If  a  flame  is  not  established  on  the  igniter  within  10  seconds, 
close  the  igniter  safety  valve  and  determine  and  correct  the  cause  of  the  ignition  failure.  Wait  at  least  1 
minute  before  a  retry.  Halt  the  startup  sequence  after  the  second  failed  attempt  to  ignite  the  igniter. 

11.  Note  that  variations  of  the  inocedure  are  requited  for  boilers  equipped  widi  Qass  3  special 
electric  igniters. 

12.  Ensure  foat  the  igniter  is  operating  cotT«:tly  and  proceed  to  supply  fuel  to  the  main  burner. 
Failure  to  ignite  within  5  seconds  will  cause  a  master  fuel  trip.  After  a  trip,  the  startup  sequence  will  be 
repeated  from  step  6.  The  second  failure  to  lightoff  will  halt  the  startup  procedure. 

13.  Wait  for  operator  verification  that  a  stable  flame  has  been  established.  Ramp  bunKr  air  registers 
or  dampers  to  die  normal  operating  posidoiL  Engage  the  fuel  flow  control  loop  widi  seqwint 
corresponding  to  die  current  fuel  flow  level.  Be  sure  that  ignition  is  not  lost 

14.  Shut  down  the  igniters  and  verify  that  a  flame  is  maintained. 
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15.  Shut  off  the  burner  header  atmospheric  valve  if  necessary  (see  NFPA  8SB  section  S-2.1.2  item 
P)- 

16.  If  the  boiler  is  equipped  with  other  burners,  place  them  in  service  as  required  and  as  specified 
in  NFPA  8SB  section  S-2.1.2  item  q. 

17.  For  multiple  burner  systems,  ignite  the  maximum  number  of  burners  consistent  with  anticipated 
continuous  load. 

18.  Verify  that  a  stable  flame  has  been  established.  Ramp  fuel  flow  to  the  minimum  operating  level 
if  necessary.  Nonoperatirig  burners  should  be  corr^y  placed  out  of  service.  Ramp  airflow  down  to  a 
level  sudi  that  the  excess  oxygen  level  is  about  10  percent  above  the  normal  operating  target  Using  this 
setting  as  the  excess  oxygen  seqxrint,  engage  the  excess  oxygen  control  loop.  Verify  burner  operaticn 
again.  Engage  automatic  excess  oxygen  setpoint  control.  Verify  burner  operation.  Use  drum  piossure 
as  the  pressure  loq)  controlled  variable  and  close  the  steam  pressure  loop.  Ramp  up  the  drum  pressure 
to  match  the  steam  header  pressure.  Switch  the  steam  pressure  controlled  variable  to  header  pressure  to 
place  the  boiler  on  line. 

Shutdown  Procedure 

The  shutdown  procedure  complies  with  NFPA  8SB  section  S-2.3.  The  following  generalized 
procedure  will  be  customized  as  required  for  individual  installations: 

1.  Switch  the  pressure  loop  controlled  vatiidtle  to  drum  pressure  to  t»ing  the  boUer  off  line. 

2.  Disengage  the  pressure  loop. 

3.  Ramp  down  the  fuel  flow  setpoint  until  tl»  airflow  reaches  the  purge  airflow  rate. 

4.  Disengage  the  excess  oxygen  control  Io(^  and  set  the  air  flow  rate  setpoint  at  die  purge  airflow 
rate. 

5.  Continue  to  ramp  the  fuel  flow  down  to  a  level  corresponding  to  startup  conditions  and 
disengage  the  fuel  flow  loop.  The  burners  should  txrw  be  operating  at  the  startup  fuel  flow  setting.  IMace 
burner  dampers  in  the  correct  startup  position.  Shut  down  each  burner  by  closing  its  safety  shutoff  valve. 
Shutdown  die  last  taimer  by  tripping  the  main  safety  shutoff  valve.  Close  all  burner  safety  shutoff  valves. 
Open  all  atmosidieric  vent  valves. 

6.  Cmnplete  a  unit  purge. 

7.  Shut  off  the  feedwater. 

8.  Shut  down  die  blowers  and  opdotudly  close  the  dampers. 

Flow  Charts  of  the  Automatic  Startup  and  Shutdown  Procedures 

Figures  15  to  24  show  die  flow  charts  of  the  starhip  and  shutdown  procedures. 
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F^ure  15.  Flow  Chart  for  Main  Startup/Shutdown  Procedure. 


Figure  16.  Flow  Chart  of  Startup  Subroutine. 
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Figure  17.  Flow  Chart  of  Startl  Subroutine. 
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Figure  18.  Flow  Chart  of  Start2  Subrouthie. 
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Figure  19.  Flow  Chart  of  SUu13  Subroutine. 
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RETURN. 


Figure  21.  Flow  Chart  of  Starts  Subroutine. 
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END 


Figure  22.  Flow  Chart  oi  Shutdown  Subroutine. 


48 


Figure  23.  Flow  Chart  of  Shutl  Subroutine. 
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niOMPT:  lYess  PgDn  to  disengage  FF  c<»itroller 


dRETURN^ 


Figure  24.  Flow  Chart  of  Shut!  Subroutine. 
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Implement  the  Automatic  Startup  and  Shutdown  Procedures 


FIX  DMACS  Grairfiic  Bivironment.  FIX  DMACS  provides  a  graphical  environment  that  allows 
tte  user  to  monitor  and  control  industrial  processes  using  custom-designed  gr^hical  representations  ot 
the  process.  The  graphical  representation  can  be  either  a  simple  diagram  or  a  detailed,  animated  display. 
Each  display  contains  links  that  display  messages,  values,  alarms,  and  real-time  trends  of  the  process.  The 
FIX  DMACS  gra}^c  environment  developed  for  boiler  automatic  startup  and  shutdown  consists  of  the 
foUowing  screens:  BOILER  STARTUP  &  SHUTDOWN.  INNER  LOOP  PID,  OUTER  LOOP  PID,  and 
OUTPUT.  All  of  these  screens  are  linked  together  in  the  following  order;  BOILER  STARTUP  & 
SHUTDOWN  ->  INNER  LOOP  PID  ->  OUTER  LOOP  PID  -->  OUTPUT  ->  BOILER  STARTUP  & 
SHUTDOWN  (Appendix  D). 

The  BOILER  STARTUP  &  SHUTDOWN  screen  is  the  master  screen  for  boiler  automatic  startup 
and  shutdown.  In  th's  screen,  the  options  STARTUP  or  SHUTDOWN  is  given  by  a  MAIN  MODES  field 
choice,  and  a  sequence  of  prompts  necessary  for  the  startup  or  shutdown  procedure  are  given  (Appendix 
E).  Moreover,  all  the  setpoints  of  the  p:''*ess  variables,  WF  (water  flow),  AF  (air  flow),  FF  (fuel  flow), 
WL  (drum  water  level),  DP  (drum  pressure),  and  Oj  (stack  excess  oxygen),  can  be  assigned  on  the  screen. 
Furthermore,  all  the  controllers  can  be  engaged  (or  disengaged)  by  pressing  the  PgDn  or  PgUp  key. 

The  INNER  LOOP  PID  performs  PID  control  for  actuators,  i.e.,  it  controls  the  water  flow  (WF), 
air  flow  (AF),  and  fuel  flow  (FIO  loops.  In  this  screen,  all  the  controller  parameters  can  be  easily  set  from 
the  display.  Also,  all  the  seqroints  of  the  actuators  can  be  set  on  the  display  as  well.  The  OUTER  LOOP 
PID  screen  performs  the  same  function  as  the  INNER  LOOP  PID,  except  it  controls  the  outer  loop 
variables  water  level  (WL),  drum  pressure  (DP),  and  stack  excess  oxygen  (02).  The  OUTPUT  screen 
graphically  shows  all  the  seqxrints  and  measured  values  of  the  process  variables  in  both  numerical  and 
bar  gra{^  format,  and  shows  the  historical  trend  of  the  setpoints  and  measured  values  of  the  drum 
pressure,  stack  excess  oxygen,  drum  water  level,  and  steam  flow. 

Operating  Procedure.  To  perform  the  automatic  startup  and  shutdown,  the  equipment  required  must 
include  a  )iMac-6(X)0  boiler  controller  as  well  as  an  IBM  compatible  PC.  The  PC  should  have  boiler 
control  software  installed  in  a  directory  called  BOIL,  and  the  FIX  DMACTS  package  installed  in  the 
DMACS  directory.  The  following  steps  outline  the  process: 

a.  Coimect  the  COMl  port  of  the  386SX  PC  to  the  COMO  port  of  pMac-6000.  Turn  on  the 
power  of  the  pMac-6000. 

b.  Change  to  the  BOIL  subdirectory  of  the  PC.  This  directory  should  contain  the  Kermit  file 
transfer  protocol  and  the  executable  boiler  control  code.  Type  KERMIT  to  establish  the  communications 
between  the  386SX  PC  and  the  pMac-6000. 

c.  Type  "SET  BAUD  19200”  to  define  the  baud  rate.  £>ownload  the  boiler  control  executable 
file  to  the  pMac-6000  using  Kermit 

d.  Press  CTRL-]  (press  “CTtrl”  and  “]”  keys  at  the  same  time),  then  i»ess  “C”  key  to  exit 
KERMIT.  Type  “EXIT’  to  return  to  the  BOIL  subdirectOTy. 

e.  Change  to  the  DMACS  subdirectOTy.  Type  “DMACS”  to  start  the  FIX  DMACS  program. 

f.  CThoose  the  SAC4  function  in  the  background  menu.  Type  “Alt-x”  to  start  tire  scan  for  node 

BOILER. 
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g.  Qioose  the  PVIEW  function  in  foreground  menu.  Type  “NEWBLR”  to  load  FIX  DMACS 
graphic  environment  A  screen  named  BOILER  STARTUP  &  SHUTDOWN  will  appear  on  the  display. 

h.  Move  the  cursor  to  the  MAIN  MODES  field;  type  “1”  for  startup  or  “2”  for  shutdown.  A 
prompt  STARTUP  NOW!  or  SHUTDOWN  NOW!  will  ^pear  on  the  screen. 

i.  Move  the  cursor  to  the  field  PROMPT  RESPONSE,  the  program  will  display  a  sequence 
of  prompts  during  execution  of  the  automatic  startup  or  shutdown  procedure.  Type  “0”  in  the  PROMPT 
RESPONSE  field  to  move  to  the  next  step  when  the  prompted  task  has  been  completed. 


Furnace  Draft  Control  Loops 

The  draft  is  defined  as  a  “current  of  air.”  In  balanced  draft  boiler  systems,  the  furnace  pressure  is 
designed  to  be  slightly  negative  relative  to  atmospheric  pressure  to  ensure  that  any  leakage  will  result  in 
relatively  cool  combustion  air  leaking  into  the  furnace  instead  of  very  hot  combustion  gases  leaking  out 
Such  boiler  systems  usually  lely  on  the  use  of  an  induced  draft  fan  in  combination  with  a  forced  draft  fan. 
The  induced  draft  fan  is  used  to  reduce  the  furnace  pressure  and  to  ensure  that  it  is  always  negative  with 
respect  to  atmospheric  pressure.  In  normal  practice,  the  furnace  pressure  or  draft  is  controlled  to  a  very 
slightly  negative  pressure  setpoint  by  regulating  either  or  both  the  forced  and  induced  draft  fans.  In  the 
balanced  draft  boiler,  the  forced  and  induced  draft  fans  share  the  load  of  moving  the  combustion  air  and 
flue  gases  through  the  system.  The  balance  point  is  the  pressure  or  draft  in  the  furnace.  This  inessure 
level  is  determined  by  the  relative  amounts  of  “push”  and  “puli”  of  the  forced  and  induced  drafts, 
respectively.  The  ftunace  draft  controller  involves  two  control  loops:  the  windbox  pressure  control  loop 
and  furnace  pressure  control  loop.  The  object  is  to  control  the  speed  of  the  forced  and  induced  fans  such 
that  the  furnace  pressure  is  maintained  at  a  slightly  negative  setpoint  Conventional  control  methods  for 
draft  control  mainly  are  PID  control.  Figure  2S  shows  the  control  logic  diagram.  . 

Since  the  system  is  subject  to  dynamic  load  change  from  the  steam  pressure  demand,  and  the  fur¬ 
nace  draft  measurement  is  subject  to  ccmsiderable  process  “noise,”  it  is  very  difficult  to  tune  the  PID 
controller  to  get  satisfactory  results  (Dukelow  1986).  The  prqwsed  Generalized  Prediction  Control  (GPC) 
schemes  for  the  Furnace  pressure  control  loop  and  Windbox  pressure  control  loop  (Figures  26  and  27) 
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Figure  25.  Draft  Control  f  ogle  Diagram. 
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Figure  26.  Flimace  Pressure  GPC  Control  Loop  Diagram. 
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Figure  27.  Windbox  Pressure  GPC  Control  Loop  Diagram. 


should  improve  system  performance.  Computer  simulations  used  the  proposed  Generalized  Prediction 
Control  (GPC)  algorithm  to  compare  the  p^otmance  by  die  PID  ctmtrol  and  the  GPC  control  based  on 
an  artificial  model  of  the  processes.  The  model  parameters  were  determined  by  off-line  identification 
using  the  input-output  data  acquired  from  the  plant  (AstrOm  and  Wittenmaile). 

Furnace  Pressure  Control  Loop 

The  fiimace  pressure  model  assumes  that: 


A(q')y(lc)  =  q'"B(q"‘u(k)  +  C(q‘)e(k) 


Where 

y(k)  is  die  fiimace  pressure 

u<k)  is  die  controller  output,  the  induced  draft  fan  qieed  command  signal 
e(k)  is  die  dimnbance  fi^  plant  steam  pressure  demand. 


[Eq52] 


Assume  that  the  polynomials  are: 


A(q-')  =  l-.8959q-'-0.043q-' 

Bl(q-')  =  0.2778  +  1.503q-' 

C(q-')  =  1 

Let  m  =  0.0(X2,  and  the  disturbance  signal  e(k)  firom  the  steam  demand  be 

e(k)  =  0.l[sin(o(k)+sin(3o)k)-^sin(8<ok)+sin(lSa)k)+sin(SO(i)k)] 
as  shown  in  Figure  28. 

Figure  29  shows  the  furnace  pressure  output  using  the  GPC  control  scheme,  where  the  desired 
furnace  pressure  is  set  as  -0.S  in.  of  HjO.  It  is  clearly  shown  that  despite  the  presence  of  the  large 
disturbance  (the  signal-to-noise  ratio  is  about  0.5),  the  output  of  the  control  system  is  close  to  the  setpoint. 

When  a  PID  control  scheme  is  used,  the  system  performance  is  greatly  influenced  by  the 
disturbance.  If  there  is  no  disturbance,  the  PID  control  system  can  follow  the  seqwint  very  well  (Figure 
30).  A  disturbance  however,  will  deteriorate  system  performance.  Hgure  31  shows  the  furnace  pressure 
output  with  PID  c(Hitrol,  where  the  disturbance  level  is  the  same  as  with  GPC  control.  The  furnace 
pressure  ouqiut  oscillates  with  large  tracking  errors  due  to  the  large  disturbance. 

Windbox  Pressure  Control  Loop 

The  simulation  model  for  windbox  pressure  loop  assumes  that: 

A(q-')y(k)  =  q^B,(q-')u,(k)  +q-^B,(q-')u2(k)  +C(q-')e(k)  fEq  53] 


Where 

y(k)  is  the  windbox  pressure 

u,(k)  is  the  controller  output,  which  is  the  fenced  draft  fan  speed  command  signal 

Ujfk)  is  the  feed  forward  signal  of  Oj 

e(k)  is  the  disturbance  from  the  plant  steam  pressure  demand. 

The  polynomials  are: 

A(q-')  =  l-0.8959q-'-.043q-^ 

B,(q->)=.2778+1.503q-' 

B2(q-')  =  1 
C(q-')  =  1 


and  loop  delay  n,=l  second,  njsO.  Assume  that  die  feedforward  signal  of  O2  is  6  percent,  and  the  desired 
setpoint  of  windbox  pressure  is  1  in.  of  HjO.  Figure  28  shows  the  disturbance  level  from  the  steam 
pressure  demand.  Figure  32  shows  the  windbox  pressure  output  when  a  PID  controller  was  used.  Because 
of  the  large  disturbance  (the  ratio  of  the  signal  to  noise  is  0.5),  the  POD  control  cannot  keep  tracking  the 
seqioint,  so  the  system  performance  is  unacceptable.  Figure  33  shows  the  result  with  the  GPC  controller. 
De^te  die  presence  of  the  large  disturbance,  the  windbox  pressure  is  maintained  close  to  the  setpoint  1 
in.  td  HjO. 
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Figure  31.  Furnace  Pressure  Output  With  PID  Control. 
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Figure  32.  Windbox  Pressure  Output  With  PID  Control. 


Figure  33.  Windbox  Pressure  Output  With  GPC  Control. 


4  ABBOTT  PLANT  TEST  RESULTS  AND  DISCUSSIONS 


System  testing  was  peiformed  in  the  summer  of  1991.  The  goal  of  the  test  program  was  to  demon¬ 
strate  that  the  prototype  control  system  was  able  to  successfully  cxmtrol  an  operational  boiler.  The  test 
results  wcHild  provi^  data  for  comparison  of  performance  of  the  prototype  controller  with  conventional 
control  systems. 

The  subsequent  tasks  involve  incmporating  several  enhancements  to  die  control  system  to  bring  it 
closer  to  the  level  of  performance  needed  for  extorsive  field  trials.  These  tasks  include  enhancements  to 
the  software  to  permit  automatic  startup  and  shutdown  of  the  controlled  boiler,  and  addition  of 
“robustificaticHi’’  schemes  (Bitmead,  Gevers,  and  Wertz  1990;  Middleman  and  Goodman  1990)  to  the 
model  process  identification  algorithm  used  by  the  Generalized  Predictive  Controller  (Qarke  1988;  Clarke, 
Mtdita^,  and  Tuffs  1987).  These  robustification  schemes  are  designed  to  ensure  that  the  controller  will 
cerate  well  even  under  seldmn  encountered  condidcsis.  Finally,  certain  schemes  designed  to  permit  safe 
boiler  operatimi  in  case  of  degraded  boiler  conditions  as  well  as  techniques  for  diagnosis  of  boiler 
(grating  equipment  failure  or  degradation  were  U>  be  identified. 


Identification  of  System  Defidendes  and  Remedies 

The  current  control  system,  based  on  an  Analog  Devices  pMac  6000  (Analog  Devices  1988, 1991) 
controller  with  an  IBM  PC  type  286  or  386  host  computer  miming  the  FIX  software  package,  continues 
to  serve  as  a  successful  development  vehicle,  even  though  several  limitations  have  become  apparent 
Many  of  these  limitations  are  a  direct  result  of  the  fact  that  state-of-the-art  commercial  control  equipment 
has  not  been  designed  to  accommodate  requirements  of  leading  edge  control  schemes  such  as  the 
Generalized  Predictive  Omtrol  algorithm.  The  major  limitations  are: 

1.  Communications  capabilities  between  the  PC  running  the  FIX  as  the  operator  interface,  and  the 
pMac  running  the  control  co^  are  limited.  The  amount  of  data  that  can  be  passed  as  one  block  of 
information  is  limited  to  100  items.  Of  more  concern  is  the  fact  that  variables  must  be  prescaled  into  a 
fixed  range.  This  feature  is  designed  to  accommodate  the  needs  of  fixed  range  transducers  and  actuators, 
but  it  has  been  a  major  limitation  when  trying  to  transfer  data  such  as  die  parameters  in  the  plant  model 
used  by  die  GPC  algorithm. 

2.  The  moderate  computing  power  of  the  pMac  60(X)  is  a  limitation.  As  reported  elsewhere  in  this 
document,  a  higher  sampling  rate  will  be  necessary  to  correcdy  control  boilers  of  die  ctqiacity  in  use  at 
Abbott  Power  Plant  In  addition,  several  enhancements  to  the  GPC  controllers  have  bem  added.  Bodi 
factos  will  place  a  heavier  computing  load  cm  the  pMac  6(X)0.  It  is  not  clear  yet  if  the  pMac  can  meet 
diese  computing  demands.  It  is  clear  that  diere  is  a  limit  that  will  most  probably  be  reached  at  some 
point  It  must  be  emjdiasized  that  as  yet  no  computing  limitatimis  have  been  encountered.  The  pMac  can 
easily  cmitrol  die  boiler  when  configured  as  a  conventional  PID  controller. 

3.  The  pMac  6(X)0  has  two  unfortunate  functimis  diat  occur  during  a  software  upset  (H*  reset.  The 
most  important  is  the  switching  of  all  control  outputs  to  minimum  output  level.  Normally,  outputs  diould 
be  held  at  their  last  valid  setting  after  a  software  upset.  The  second  is  the  inability  of  the  pMac  to  execute 
an  orderly  restart  after  an  upset  Due  to  these  limitations,  the  pMac  cannot  by  itself,  be  entrusted  with 
boiler  combustion  control. 

4.  The  pMac  6(XX)  does  nm  have  a  built-in  power  failure  prevention  system. 
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S.  The  memoty  used  by  ttie  |iMac  6000  lacks  error  detection  and  correction  circuitry.  This  means 
diat  neither  data  nor  code  stored  in  random  access  memory  (RAM)  is  (Hrotected  against  caused  by 
bad  memory  locations,  <x  bits  lost  due  to  cosmic  radiation  and  power  fluctuations. 

Some  limitations  described  above  have  already  been  addressed.  Solutions  to  each  limitation  will 
be  discussed  below. 

1.  The  limited  communications  ability  of  the  FIX  software  package  will  be  addressed  by  switching 
to  Intelhition’s  newer  software  package,  the  FIX  DMACS  (Intellution  1991).  This  package  makes  use  of 
much  better  ccunmunication  software  and  incorporates  many  other  useful  enhancements.  The  switdi  from 
die  FIX  to  die  FIX  DMACS  is  straightforward.  DMACS  is  delivered  with  tools  to  help  convert  a|q>lica- 
tions  from  die  FIX  to  the  FIX  DMACS. 

2.  The  moderate  computing  power  of  the  pMac  6000  is  a  continuing  concern.  The  stated  goal  of 
diis  fnoject  has  been  to  produce  a  set  of  control  algoridims  and  mediods  that  could  be  easily  ptxted  to 
odier  platftmns.  As  a  result,  the  computing  limitations  of  the  pMac  6000  are  not  seen  as  a  major 
roadblock.  If  die  conqiuting  crqiacity  of  die  pMac  6000  is  reached,  anodier  platform  will  be  us^ 
Therefore,  a  search  ftv  an  alternate  control  device  has  begun.  At  inesent,  the  most  pitunising  candidate 
is  GE  Fanuc’s  Series  90-70  Model  781  PLC.  This  device  is  a  high  quality  programmable  logic  controllo’ 
widi  optimal  C  programming  capability.  The  process^  on  this  device  is  an  Intel  80386  tunning  at  twice 
the  clock  speed  of  the  Intel  80188  on  the  pMac  6(XX).  Clearly,  this  device  should  (xovide  a  large  boost 
in  computing  power.  In  addition,  many  interlocks  required  by  NFPA  Standards  8SA,  8SB,  and  8SG  can 
be  implement^  using  the  ladder  logic  capability  of  the  90-70.  The  90-70  uses  a  special  purpose  custom 
processor  to  execute  die  ladder  logic,  so  this  function  would  have  limited  impact  on  the  cmtrol  algmithms. 
The  soies  90-70  is  apparendy  die  «ily  currently  available  high  quality  I^  which  includes  a  general 
purpose  computing  (^on  adequate  to  meet  the  requiremmts  of  our  system. 

3.  For  testing  purposes.  Control  Technology,  Inc.  Model  7312  Lotqmiate  protection  devices  have 
been  placed  between  the  pMac  controller  and  die  boiler  fuel,  air,  and  water  flow  actuators.  Upon  a  pMac 
watchdog  timer  time-out,  the  actuator  command  signals  are  frozen  at  the  values  they  had  S(X)  ms  before 
the  time-out  These  devices  also  permit  manual  control  of  die  boiler  from  their  frcnt  panel. 

The  Loopmates  (xobably  could  serve  as  a  long-term  solution  to  die  problem  caused  by  Mac  outputs 
switching  to  low  during  an  upset  A  mme  serious  problem  is  posed  by  NFPA  Standard  85B,  “Prevention 
of  Furnace  Explosions  in  Natural  Gas-Fired  Multiple  Burner  Boiler-Furnaces  1989  Edition,"  section  4- 
2.2.3,  which  states  that  '‘The  logic  system  performing  the  safety  functions  for  burner  management  shall 
not  be  combined  widi  any  odier  logic  system.”  This  means  that  a  second  control  device  must  be  dedicated 
to  bumo’  management  This  standard  may  not  apply  to  some  central  heating  plants.  NFPA  8SA, 
“Prevention  Furnace  Explosions  in  Fuel  Oil-  and  Natural  Gas-Fired  Single  Burner  Boiler-Furnaces, 
1987  Edition”  may  rqiply  to  many  central  heating  plants.  NFPA  8SA  does  not  require  an  indqiendent 
burner  management  controller.  The  requirement  of  NFPA  8SB  is  prudent  and  because  a  standardized 
control  system  is  to  be  designed,  it  would  appear  wise  to  satisfy  die  requirements  of  NFPA  8SB. 

As  a  result,  it  may  be  best  to  insert  a  progiaminable  logic  ccmtroller  (PLC)  between  the  pMac  and 
the  boilm'.  The  PLC  u^I  monitor  and  pass  through  the  fuel  and  air  command  signals.  As  required  by 
section  4-2.2.3  of  NFPA  85B,  die  PLC  will  be  i»t>grammed  to  implement  purge  interlocks  and  timing, 
mandatcny  safety  shutdowns,  trial  timing  fOT  ignition,  and  flame  monitoring.  In  additicm,  die  device  wUl 
be  pfogtaummed  to  test  for  dangerous  air^l  ratios  and  to  lock  die  air  and  fuel  command  signals  upon 
a  )Afac  watchdog  timer  upset  The  PLC  will  itself  be  chosen  to  meet  the  requiremoits  of  NFPA  8SB, 
including  die  requirements  of  sections  4-2.2.1  and  4.2.2.2.  At  present,  it  would  aj^pear  that  a  simple  PLC 
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can  meet  tfie  requirements  of  NFPA  8SB.  One  good  candidate  is  the  GE  Fanuc  Series  90-30  Model  31 1 
PLC. 


GE  Fanuc  has  fiuther  expanded  its  Series  90-30  PLC  CPU  family  with  the  addition  of  the  mid-range 
Model  341.  The  new  341  CPU  includes  an  80188-based  processor  with  special  ct^rocessor  that  operates 
at  an  accelerated  clock  rate  of  20  MHz  compared  to  the  previous  8  MHz.  The  90-30’s  memory  has  been 
enlarged  by  a  facttn*  of  five,  to  40K  words.  GE  Fanuc  also  has  increased  the  number  of  registers  from 
2,000  to  10,000.  The  increased  clock  speed  boosts  the  Boolean  I/O  scan  rate  from  0.4  msec/K  to  0.3 
msec/K.  And  non-Boolean  operations  occur  in  half  the  time  previously  required.  Analog  inputs  for  the 
90-30  have  increased  fiom  128  to  1024,  while  its  analog  outputs  have  climbed  from  64  to  2S6.  This 
increased  speed  and  IA3  capability  means  that  the  90-30  can  be  used  in  high-speed  distributed  sqiplications 
such  as  high  speed/Iarge  quantity  material  handling,  which  previously  required  a  more  expensive  and 
larger  H.C. 

GE  Fanuc’s  high-end,  open-architecture,  90-70  PLC  also  improved  its  CPU  performance.  One  big 
advantage  to  using  open  systems  and  irKiustry  starrdards  is  the  relative  ease  of  upgrading  system  processors 
and  support  peripherals  while  maintaining  a  seamless  upgrade  path. 

The  90-70  now  includes  a  80486DX-based  Model  914  CPU.  As  a  result,  this  CPU  has  the  same 
Boolean  scan  rate  (0.4  msec/K)  as  other  90-70  CPUs,  but  now  includes  the  capability  to  handle  higher 
level  math  functions  more  quickly,  and  features  an  expanded  user  memory  of  S12K  bytes  (GE  Fanuc 
Automation  1993).  To  help  ease  system  maintenance,  the  Model  914’s  operating  system  can  be 
downloaded  from  disk  to  flash  PROM,  to  ease  and  simplify  O/S  updates.  The  Model  914  is  intended  for 
more  complex  control  applications  and  can  be  {»ogramm^  in  relay  ladder  logic,  state  logic,  sequential 
function  block  (SFC),  and  in  C — like  the  other  90-70  CUPs — to  suit  individual  program  development 
needs. 

The  Model  914  CPU  is  capable  of  industry-standard  conununicadons  protocols,  including  Ethernet, 
MAP,  RTU/Modbus,  CCM,  SNP  and  GE  Fanuc’s  Genius  intelligent,  distributed  VO  system,  and  the  I/O 
link.  The  last  provides  an  interface  between  GE  Fanuc  PLCs  and  its  computer  numerical  ctmtrols.  Even 
if  the  pMac  is  replaced  with  a  programmable  logic  controller  with  computational  capability  (see  above), 
a  second  burner  management  controller  will  be  required. 

4.  The  fact  that  the  pMac  6000  (and  the  Loopmate  devices  currently  used  to  ensure  that  the 
actuatin’  conunand  signals  are  held  at  reasonable  values  after  an  upset)  does  not  have  inherent  power 
failure  protection,  should  not  be  a  cause  for  much  concern.  The  simple  addition  of  an  uninteimptable 
power  supply  (for  example,  a  Superior  Electric  Company  Model  UPS6100SR  Stabiline)  can  easily  ensure 
control  system  operation  for  the  required  time. 

5.  The  lack  of  memory  error  detection  and  correctitm  circuitry  on  the  pMac  controller  is  of  scnne 
concern.  If  the  pMac  is  replaced  with  a  more  powerful  computing  device,  one  of  the  selection  criteria  for 
the  replacemoit  device  will  be  memory  protection.  The  problem  is  not  as  severe  as  it  might  first  appear. 
In  die  first  place,  the  addition  of  a  line  cmiditicmer  and  uninterruptable  power  supply  such  as  die  unit 
discussed  above  will  greatly  lessen  the  likelihood  of  memory  corruption  due  to  power  fluctuations.  While 
the  RAM  memory  of  the  pMac  has  neither  error  correctitm  nor  error  detection  circuitry,  it  is  static  and 
not  dynamic  RAM.  Static  RAM  is  far  more  resistant  to  corruption  due  to  power  fluctuation  arul  high 
energy  particles  than  is  dynamic  RAM.  Finally,  the  addition  of  a  burner  management  controller  as 
required  by  NFPA  85B  along  with  the  pMac’s  watchdog  timer  crqiability  provides  a  measure  of 
redundancy  and  protection  against  failure  of  the  pMac  controller  due  to  memory  corruption. 
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Tcsiilig  of  tile  Prototype  Control  System 
Outline  cf  the  Test  Plan 


The  following  is  an  outline  of  the  test  sequence  at  the  Abbott  Power  Plant  boiler  no.  2: 

1.  Stabilize  the  boiler  at  an  acceptable  operating  point.  The  excess  oxygen  level  would  be  set  high 
as  a  safety  measure. 

2.  The  fuel,  air,  and  water  flow  actuator  loops  were  tuned  next.  These  lo<q)s  are  ccmventional  PID 
control  loops  implemented  digitally  on  the  Mac  controller. 

3.  Once  the  actuator  control  loops  were  tuned,  each  outer  loop  control  model  was  identified,  and 
after  identification,  the  controller  was  started.  The  fust  outer  loop  model  to  be  identified  was  the  water 
level  loc^.  Afto*  ctmvergence  of  this  model,  the  GPC  controller  fcM*  water  level  was  started.  Next,  this 
procedure  was  repeated  for  excess  oxygen.  Finally,  the  drum  pressure  control  model  and  controller  was 
started. 

4.  Next  small  perturbations  of  setpoint  for  each  control  loop  were  af^lied  to  test  transient  response. 

5.  A  steam  deiiuuid  disturbance  was  simulated  by  changing  the  fuel  rate  on  boiler  no.  3  (on  the 
same  steam  header). 

6.  A  large  change  in  steam  demand  was  initiated  next  by  making  a  large,  gradual  change  in  the 
fuel  rate  of  boiler  no.  3.  As  a  safety  measure,  the  GPC  controller  was  inactive  during  the  transition.  The 
boiler  was  in  manual  control  mode  during  the  transition.  After  reaching  the  new  operating  point,  steps 
2  through  S  were  repeated.  The  {vinciple  goal  of  this  procedure  was  to  determine  the  chwges  in  the 
controlier  riKxlel  at  different  operating  points. 

7.  The  boiler  was  tetumed  to  the  original  set  of  operating  conditions  widi  the  controllers  on.  This 
experiment  demonstrated  die  self-tuning  tracking  ability  of  the  controllers. 

8.  The  perturbation  experiment  was  repeated  at  this  operating  point. 

9.  Finally,  this  entire  test  plan  was  repeated  using  digital  PID  controllers  to  create  a  comparison 
data  set 

Appendix  A  gives  details  of  dus  test  plan. 

Test  Difficulties 

Problems  Tuning  Actuator  PID  Loops.  The  parameters  of  the  PID  controller  were  tuned  by  using 
standard  tuning  rules  from  Astrdm  (1984)  (Appradix  B).  Using  the  transient-response  method,  the 
contndlers  worked  as  expected  except  that  the  controlled  variables  contained  limited  cycle  oscillations. 
The  pneumatic  PID  controllers  nomially  used  in  Abbott  also  produce  such  oscillations  in  the  ccmtrolled 
variable,  but  at  a  smaller  amplitude.  Much  time  was  spent  in  an  effort  to  satisfy  die  operators  that  actuator 
control  was  acceptable.  The  solution  to  the  larger  dian  normal  limit  cycles  is  probably  a  faster  sampling 
rate. 


Problems  With  the  Identifier  on  the  Water  Level  Loon.  It  was  found  to  be  impossible  to  obtain 
convergence  for  plant  model-parametos  associated  widi  die  water  level  loop.  Since  the  GPC  is  an 
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adaptive  mcxtel-based  method,  the  effoits  to  use  the  GPC  cm  the  water  level  loop  were  unsuccessful. 
Instead,  a  fixed  plant  model  with  predictive  control  for  the  water  level  loop  was  used  for  the  test. 

Test  Results 

Drum  Pressure  Loop.  The  drum  pressure  loop  performed  well  under  GPC  control  (Figure  34). 
Though  the  setpoints  were  varied  over  a  wide  range  (from  323  to  347  psi  in  less  dian  1000  seconds),  the 
loop  output  tracked  die  setpoint  nicely.  This  is  quite  difficult  for  conventional  pneumatic  PID  control. 
Tij^ter  drum  pressure  control  can  improve  upset  recovery  and  operational  safety.  The  Abbott  Power  Plant 
personnel  all  agreed  that  the  drum  pressure  control  was  superior  using  the  GPC  algorithm. 

Excess  Oxygen  Loop.  In  general,  the  perfcHmance  of  the  excess  oxygen  loop  under  GPC  control 
is  excellent  (Figure  35,  data  batch  2000-6000  seconds).  During  the  sharp  reduction  in  boiler  steam  flow 
diat  occurred  after  6000  seccMids,  control  degraded  ccHisiderably.  It  was  found  that  this  prc^lem  mainly 
was  caused  by  a  deficiency  of  die  identifier.  The  excess  oxygen  loop  is  nonlinear  so  the  model  dynamics 
change  at  different  operating  regimes.  To  check  this,  <me  mi^t  choose  a  batch  of  data  from  5001  to  7500 
seconds,  and  use  the  PEM  batch  parameter  estimation  function  in  MATLAB  to  estimate  the  model 
parameters.  Under  these  conditions,  the  maximum  cme  step  ahead  predictitMi  error  is  about  2.8  percent 
of  oxygen  (Figure  36).  By  shifting  the  data  set  from  6(X)1  to  8500  seconds,  however,  the  maximum  one- 
step-ahead  prediction  error  would  be  only  about  1.6  percent  of  oxygen  (Figure  37).  CHearly,  the  model 
parameters  needed  to  change  sharply  after  6000  seccmds  into  the  test  Since  a  slow  sampling  rate  (5 
seconds  fm*  excess  oxygen  loop)  was  used  during  dte  test  the  identifier  could  not  quickly  catch  the  rapid 
changes  of  the  parameters  in  the  model  so  that  die  predictioi  control  algorithm  was  no  longer  able  to 
control  the  excess  oxygen  level  satisfactorily.  Using  the  data  set  logged  during  the  test  it  has  been 
discovoed  that  the  identifier  will  track  the  parameters  well  when  the  process  is  sampled  at  a  1-second 
update  rate. 

Water  Level  Loop.  Test  results  fix’  this  loop  are  marginal  (Figure  38),  since  a  fixed  model  and  a 
predictive  controller  were  used.  Nevertheless,  the  performance  of  the  GPC  cmitroller  is  at  least  as  good 
as  the  ctxitinuous  pneumatic  PID  controller  in  the  Abbott  power  plant  Rgure  39  shows  the  results  using 
PID  controUCT.  ^ffer  the  test  analysis  shows  that  parameter  estimates  will  converge  after  40  minutes  at 
a  sampling  rate  at  20  seconds  (Figure  40).  At  a  faster  sampling  rate,  the  parameter  estimates  converge 
much  quicker.  Figure  41  shows  that,  at  a  1-second  sampling  rate,  the  estimated  parameters  converge  after 
12  minutes. 

Oscillatory  Behavior  of  Actuators.  The  fuel,  air,  and  water  flow  actuator  loops  demtmstrate 
oscillatory  behavior  (Figure  42).  The  main  cause  for  this  was  probably  the  backlash  in  the  drive 
mechanism.  The  oscillatimi  was  also  present  when  using  the  Abbott  Power  Plant  ctmtinuous  pneumatic 
HD  controllers,  but  at  a  smaller  amplitude.  A  fast  update  rate  would  reduce  the  amplitude  of  this 
oscillation,  but  it  is  likely  that  only  inqxoved  actuation  equipment  could  eliminate  such  a  hard 
nonlinearity. 

Change  in  Time  Delay  With  Operating  PoinL  Time  delay  is  an  important  factor  in  choice  of  model 
structure.  One  effective  approach  to  determining  the  process  time  delay  is  to  observe  the  step  response 
of  die  ixocess.  In  the  case  of  die  excess  oxygen  loop,  it  seems  that  die  process  delay  is  slightly  changed 
at  different  operating  points.  The  change  was  surprisingly  minor  (only  1  to  2  seconds);  a  large  time  delay 
change  had  beoi  anticipated.  It  tqipears  that  gas  buoyancy  dominates  the  gas  flow  rate  in  the  boiler.  For 
die  data  batch  6501-70(X)  secmi^,  die  correct  time  delay  was  found  to  be  13  seconds  for  fuel  flow 
seqwint  to  oxygen,  and  8  seconds  for  air  flow  setpoint  to  oxygen  (Figure  43).  If  the  time  delay  is 
inoeased  to  14  seconds  for  fuel  flow  setpoint  to  oxygen.  Figure  44  shows  that  the  shape  of  the  step 
response  becomes  unreasonable.  Further,  for  data  batch  7(X)1-7500  seconds,  the  correct  time  delay  is 
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Time  bom  Start  of  Tetl  (Seconds) 


Figure  35.  Exhaust  Oxygen  and  Setpoint  With  GPC  Control. 


Excess  Oxygen  -  %  «•  Excess  Oxygen  • 


ure  36.  One  Step  Predktion  Error  for  the  Data  Batch  5001  to  7500  Seconds. 


Figure  37.  One  Step  Prediction  Error  for  the  Data  Batch  6001  to  8500  Seconds. 
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Time  from  Stan  of  Test  (Second*) 


Figure  39.  Drum  Water  Level  and  Setpoint  With  PID  Control. 


Water  Level  - 


Figare  40.  Estimated  Parameters  and  Prediction  Error  With  20-Second  Sampling  Period. 
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Water  Level  - 


Esdmated  Parameters  a  1,  a2,  &  a3.Sampline  Isec. 
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Air  Flow  -  %  Fuel  Flow  -  %  Water  Flow  • 


Figure  42.  Oscillatory  Behavior  of  the  Actuators. 
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Figure  43.  Step  Response  of  Data  Batch  <>501  to  7000  Seconds  of  02/FF;  Delaysi3  Seconds. 


Tune  from  Suit  of  Test  (Seconds) 

Figure  44.  Step  Response  of  DaU  Batch  7001  to  7500  Seconds  of  02/FF;  Ddaysld  Seconds. 
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changed  to  14  seconds  for  fuel  flow  se^wint  to  oxygen  (Figure  4S).  If  die  chosen  time  delay  is  IS 
seconds,  the  step  lespcmse  still  looks  reasonable  (Figure  46).  If  the  time  delay  is  increased  fiutha-.  Figure 
47  shows  that  the  step  response  of  the  process  diveiges. 

Performance  Comparison  of  GPC  and  PIP.  Direcdy  comparing  the  performance  of  the  GPC 
cmtroUer  with  PID  controllers  proved  to  be  difficulL  Mainly  diis  is  due  to  the  plant  not  being 
instrumented  for  se^int  measurement  when  the  conventional  pneumatic  PID  controllers  were  in  use.  It 
was  impossible  to  match  die  operating  conditions,  but  a  limited  comparison  could  still  be  made.  Figure  48 
shows  die  exhaust  oxygen  and  setpoints  with  the  |uieumatic  PID  controller  and  the  GPC  cmitroller.  It  is 
clear  that  with  the  Grc  control,  die  step  response  of  the  exhaust  oxygen  (within  0.2  percent  Oj )  is  much 
titter  than  with  PID  control  (within  0.5  percent  Oj).  The  GPC  algcnithm  provides  excellent  cmitrol  of 
combustimi.  This,  potentially,  has  large  payback.  Ti^t  oxygen  control  reduces  losses  due  to  excess  hot 
gases  leaving  the  stack  to  a  very  low  level.  Exhaust  gas  CO  cmitrol  can  cany  this  benefit  a  step  further. 
The  system  under  development  will  be  cost  competitive  with  simple  oxygen  trim  equipment,  but  will 
permit  the  tightest  possible  combustion  control  with  the  resulting  economic  benefits. 


Figure  45.  Step  Reqionse  of  Data  Batch  7001  to  7500  Seconds  of  02/FF;  Delay=14  Seconds. 
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Figure  46.  Step  Response  of  Data  Batch  7001  to  7500  Seconds  of  02/FF;  DdayslS  Seconds. 


Figure  47.  Step  ResptHise  of  Data  Batch  7001  to  7500  Seconds  ot  02/FF;  Dday=18  Seconds. 
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Excess  Oxygen  and  Setpoint  with  PE) 


5  BOILER  COMBUSTION  EFFICIENCY  CALCULATIONS 


Combustion  efficieiicy  measurement  is  a  very  important  issue,  because  it  is  an  index  directly 
indictaive  of  boiler  performance.  It  will  be  very  helpful  for  the  operator  and/or  control  strategy  in  making 
adjustments  to  maintain  optinud  boiler  productivity.  Boiler  combustion  efficiency  is  defined  as  die  ratio 
of  the  heat  transferred  to  the  water/steam  side  cmnpaied  to  the  higher  heating  value  of  the  fuel.  The 
American  Society  of  Mechanical  Engineers  (ASI^)  has  established  the  most  widely  accepted  two 
mediods  for  bmler  efficiency  testing,  the  inputloutput  method  and  the  heat  loss  method,  which  are  obtained 
tiirou^  a  standardized  calculation  procedure  as  described  in  the  Amnican  Society  of  Mechanical 
Engineers  (ASME)  Perfmmance  Test  Code  (PiC)  4.1.  (1%S).  It  is  generally  agreed  that  the  heat  loss 
metiKxl  provides  noore  accurate  results,  while  the  input/output  method  is  more  easily  understood  and,  for 
that  reason,  is  often  preferred. 

The  basic  calculation  for  the  input/output  method  is: 


Efficiency  « 


Useful  Heat  Extracted 
Total  Heat  Input 


X  100% 


[Eq54] 


Where 

Useful  Heat  Input  =  heat  content  of  steam  to  die  header  plus  the  heat  contoit  of  blowdown  flow 
miiuis  the  heat  content  of  boiler  feedwater 

Total  Heat  Input  s  heat  released  by  the  cmnbustion. 

The  input/output  ASME  test  is  well  suited  for  gas-  and  oil-fired  boilers  where  the  fuel  firing  rates 
are  oftmi  accurately  measured  and  the  heating  values  of  the  fuels  are  known  constants.  It  is  advantageous 
to  use  because  efficiency  is  calculated  from  fuel  and  steaiT)  flowmeters.  This  medKxl  often  meshes  well 
widi  existing  plant  instrumentation,  while  disadvamages  are  that:  (1)  the  method  is  sensitive  to  flow 
measurement  errors,  and  (2)  all  fuels  must  be  accurar/'iv  -leasured. 

The  Heat  Loss  MedKxl  is  used  to  calculate  individual  losses,  totalize  them,  and  determine  boiler 
efficiency  by  subtracting  the  total  losses  firom  100  percent.  The  formula  is: 


Efficiency  =  (1  - 


)viniw. 

Total  Heat  Losses 


[Eq551 


This  equation  can  be  obtained  from  the  input/ouqxit  method  simply  by  making  the  substitution:  Total  Heat 
hiput  =  Useful  Heat  Extracted  Plus  Heat  Losses.  The  heat  loss  method  calculates  boiler  performance 
using  an  equatitm  that  is  comparatively  less  sensitive  to  measurement  errors  in  liiel  flow  or  steam  flow 
than  dun  used  in  die  input/output  mediod.  This  is  an  advantage  in  solid  fuel  firing  systems,  where  it  is 
inqiossiUe  to  obtain  a  highly  accurate  fuel  measurement  In  the  case  of  a  boiler  firing  multiple  fuels,  a 
ratio  of  die  fuel  firing  rates  is  required.  Also,  fewer  on-line  measurements  are  required  and  it  is  thought 
to  usually  provide  a  more  accurate  result  Its  main  disadvantage  is  that  not  all  losses  are  computable  such 
as  losses  to  radiation,  which  are  estitiuued  using  ASME  published  tables. 

The  foUowir^  six  heat  losses  are  those  specified  by  the  ASME  PTC  4.1  (ASME  1%S)  abbreviated 
method  in  accordance  widi  die  heat  loss  meth^  to  calculate  the  heat  losses. 
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1.  Dry  Gas  Loss.  This  loss  represents  heat  lost  in  dry  flue  gas.  The  amount  of  this  loss  depends 
on:  (a)  the  stack  gas  temperature;  the  amount  of  excess  air  used  to  Are  the  fuel;  and  (c)  the 
composition  of  the  flue  gas. 

2.  Loss  Due  to  Moisture  fron  Burning.  Hydrogen  in  hydrocarbon  fuels  forms  water  when  burned. 
The  heat  viqxMization  of  diis  water  is  lost  out  of  die  stack. 

3.  Loss  Due  to  Moisture  in  Fuel.  If  the  fuel  being  burned  contains  any  moisture,  this  wiuer  must 
be  evaporated.  The  result  is  a  heat  loss. 

4.  Loss  Due  to  Unbumed  Gimbustibles.  Unbumed  or  partially  burned  fuel  represents  a  loss. 

5.  Losses  from  Radiation  and  Convection.  These  losses  through  the  boiler  structure  are  virtually 
impossibte  to  measure.  Instead,  refer  to  ASME  published  tables  to  determine  heat  loss. 

6.  Unmeasured  Losses.  Many  minor  heat  losses  are  lumped  into  this  categmy,  using  an 
“engineering  estimate.” 

The  calculation  procedure,  eidier  input/output  method  or  heat  loss  method,  is  rathor  complicated  due 
to  the  number  of  process  variables  that  had  to  be  measured  and  the  number  of  calculations  to  be  executed. 
For  real-time  application,  it  is  better  to  find  a  simple  and  easy  way  to  perform  the  boiler  efficiency 
calculations. 

Under  the  assumption  diat  there  was  complete  combustion  and  no  water-vt^  in  the  combustion 
air,  it  is  possible  to  determine  die  combustion  efficiency  as  a  function  of  flue  gas  temperature  and 
percoitage  Oj  in  die  dry  flue  gas.  Based  on  the  data  prov^ed  by  Boiler  Efficiency  Improvement  (Maples 
and  Dyer  1981),  the  combustion  efficiency  calculation  models  for  different  fuels  can  be  obtained  by  using 
a  regression  method. 


Efficiency  of  Natural  Gas-Fired  Boiler 

Ti(%)  =  -80.6545  +  16.75665X,  -  0.45138x,^ 

+  29.4500xj  - 1 .25286x/  - 1  -0.02769T 


where: 


X,  is  the  percentage  O,  in  the  dry  flue  gas 
Xj  is  the  percentage  COj  in  the  dry  flue  gas 
T  is  the  flue  gas  temperature. 

Mmeover,  die  percentage  C02in  the  diy  flue  gas  can  be  estimated  fnmi  die  percentage  Oj  in  the  dry 
flue  gas  via  die  chemical  reaction  balance  in  ^  combustion  process.  The  percent  (X>2  for  the  natural  gas- 
flred  unit  is  calculated  by: 


CX)j(%)  =  11.80121 -0.56145Oj(%) 


The  combustion  efficiency  estimadmi  error  is  the  difference  between  the  calculated  ou^ts  and  the 
table  values  provided  in  (Maples  and  Dyer  1981).  Figure  49  shows  the  combustion  efficiency  esdmatitxi 
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errors  for  tfw  natural  gas-fired  boiler,  where  the  oxygen  is  ranged  from  0  to  10  percent,  and  the  stack 
tonperature  is  varied  fnxn  170  to  630  "F.  It  is  interesting  to  note  that  the  oxygen  value  is  usually  set  at 
the  region  from  4  to  7  percent.  In  practice,  for  this  region,  the  maximum  estimation  error  is  only  about 
0.S  percent;  even  the  temperature  varies  from  170  to  630  *’F.  Figure  SO  ^ows  a  comparison  of  the 
calculated  percent  CX>2  and  the  table  values  provided  in  Maples  and  Dyer.  It  appears  that  the  e^imated 
percent  COj  values  and  the  table  values  are  very  close. 


Efficiency  for  No.  2  Oil-Fired  Boiler 

Ti(%)  =  1202.0821 -9.651 12x,  +  0.173219x,^ 

-3.2171Xj  ♦  0.135082  +  0.590321x,Xj-0.02684T 

where  x„  x,,  and  T  are  defined  similarly  as  in  die  case  of  Natural  Gas.  The  percent  COj  for  the  no.  2  Oil 
is: 


COjCfli.)  *  15.61069 -0.7425Oj(%) 


Figures  51  and  52  show  the  combustion  efficiency  estimation  errors  for  a  no.  2  oil-fired  boiler  and  the 
comparison  of  die  estimated  percent  COj  and  the  percent  COj  table  values  from  Maples  and  Dyer  (1981). 


Efficiency  for  a  No.  6  Oil-Fired  Boiler 

T)(%)  =  23.29147 -18.3246X,  ♦  1.001373x,^ 

+  14.6999x2-0.61834x2^  0.820384x,X2 -0.027 19T 


where  x„  X2,  and  T  are  defuied  similarly  as  in  the  case  of  natural  gas.  The  percent  CO2  for  the  no.  6  oil- 
fired  unit  is  obtained  by: 

C02(%)  =  16.49536  -0.7847602(%) 


Figures  53  and  54  show  the  combustion  efficiency  estimation  errors  for  a  no.  6  oil-fired  boiler  and  the 
comparison  of  the  estimated  percent  COj  and  the  percent  COj  table  values  from  Maples  and  Dyer. 

With  the  above  equations,  the  boiler  combustion  efficiencies  can  be  easily  calculated  for  die  natural 
gas,  no.  2  oil  and  no.  6  oil  fired  boilers  in  teal  time.  Figure  55  diagrams  how  the  combustion  efficiency 
is  calculated. 

The  above  combustion  efficiency  calculation  algorithm  has  been  implemented  in  the  PC  using  the 
FIX  software  package.  The  dattdiase  used  to  calculate  the  boiler  efficiency  is  listed  in  Appendix  F. 
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F^[ure  49.  Cmnbustion  Efficiency  Estimation  Errors  for  Natural  Gas>Fired  Boiler. 
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Figure  50.  Estimated  CO,  Values  for  Natural  Gas-Fired  Boiler. 
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F^ure  51.  Combustion  Efficiency  Estimation  Errors  for  #2  Oil-Fired  Boiler. 
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Figure  52.  Estimated  Percent  CO,  for  #2  Oil-Fired  Boiler. 
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Figure  53.  Combustion  Efficiency  Elsthnation  Errors  for  #6  Oil*Fired  Boiler. 
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Figure  54.  Estimated  Percent  CO,  for  #6  Oil-Fired  Boiler. 
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Figure  55.  The  Flow  Chart  of  the  Boiler  Efficiency  Calculations. 


6  ROBUSTIFICATION  OF  THE  IDENTIFIER 


Real  industrial  systems  can  rarely  be  modeled  to  a  high  level  of  accuracy.  In  particular,  a  linear 
black-box  type  CARMA  (Ccmtrolled  Auto  Regressive  and  Moving  Average)  or  CARIMA  (Controlled 
Auto-Regiessive  and  Integrated  Moving  Average)  model  of  low  order  (usually  <=  4)  as  is  normally 
assumed  in  self-tuning  control,  can  only  be  expected  to  explain  local  system  behavior.  Real  processes  are 
fraught  with  unstructured  disturbances  and  nonlinearities.  A  boiler  is  a  highly  nmilinear  multivariable 
system  with  complex  dynamics.  When  a  low-order  linear  model  is  fitted  to  it,  the  results  will  depend 
critically  on  the  frequency  content  of  the  input  signal.  Identification  problems  in  the  presence  of  noise 
and  unmodeled  dynamics  have  been  attacked  using  two  methods,  slow  sampling  (Rohrs  et  al.  1984),  and 
prefiltering  (Astrom  and  Wittenmark  1989;  Mohtadi  1988). 


Sampling  Procedure  and  System  Noise 

The  sampling  procedure  will  automatically  introduce  a  low-pass  Hlter  (Wahlberg  1990).  Its  break 
frequency  w^(which  is  defined  as  the  frequency  at  which  the  amplitude  is  3  db  below  the  maximum  value 
of  the  amplitude),  is  proportional  to  the  inverse  of  the  sampling  period  T,  i.e.,  w^l/T.  The  bandwidth 
of  this  low-pass  filter  tends  to  infinity  with  a  converge  rate  1/T  as  T  tends  to  zero.  Therefore,  most  of 
the  high  frequency  noise  and  effects  of  high  frequency  unmodeled  dynamics  will  be  present  in  the  sampled 
data  set  if  a  fast  sampling  rate  is  used. 

A  slow  sampling  rate  is  one  of  the  possible  solutions  to  the  identification  problem  with  unmodeled 
dynamics  and  noise,  since  the  low-pass  property  of  the  slow  sampling  procedure  will  attenuate  the  gain 
at  high  frequencies.  In  the  early  development  stages  of  the  control  ^gorithm,  difficulties  were  experienced 
with  simulations  at  fast  sampling  rates.  Models  converged  to  hi^  frequency  noise  models.  For  this 
reason,  a  slow  sampling  approach  was  used  initially.  However,  slow  sampling  rate  lead  to  other  identifica¬ 
tion  problems.  The  deficiencies  observed  in  the  July  1991  tests  appear  to  be  with  identiflcation.  Slow 
sampling  rate  seems  to  be  the  primary  problem,  because  the  estimator  could  not  get  enough  real  time 
information  from  the  plant  due  to  the  data  update  rate  being  too  low.  This  resulted  in  the  parameter 
convergence  speed  being  too  low.  In  the  case  of  the  water  level  loop,  with  a  sampling  interval  of  20 
secmids,  the  parameters  required  over  40  minutes  to  converge.  The  tracking  property  of  the  identifier  was 
poor  with  a  low  sampling  rate;  there  was  a  heavy  nonlinearity  in  the  O2  loop,  and  the  identified  model 
parameters  needed  to  be  changed  with  different  operating  points.  Poor  parameter  tracking  of  the  identifier 
lead  to  the  loss  of  the  plant  model  for  the  Oj  loop  when  the  operating  points  sharply  turned  down  (Figure 
35,  6000  to  7(XX)  seconds).  Poor  identification  leads  to  poor  control. 

It  now  appears  that  better  results  can  be  achieved  by  using  fast  sampling  with  prefiltering.  When 
a  fast  sampling  rate  is  used,  precautions  must  be  taken  to  ensure  that  the  frequency  content  of  the  input 
signal  is  concentrated  in  the  frequency  range  where  the  simplified  model  is  expected  to  fit  well.  This 
implies  that  the  signals  should  be  filtered  before  they  enter  into  the  parameter  estimator. 


Filter  Types,  Break  Frequency,  and  Filter  Rational 

To  attenuate  the  effect  of  disturbances  wherever  they  are  dominant,  and  to  suppress  the  effects  of 
unmodeled  dynamics,  a  filter  must  be  chosen.  There  are  many  types  of  filters  available  for  this  purpose. 
Qassic  filters  include  Butterworth,  Bessel,  Chebyshev,  and  Elliptic  filters,  any  of  which  can  be  easily 
designed  using  standard  filter  design  methods  (Oppenheim  and  Schafer  1975;  Parks  and  Burrus  1985). 
Also,  there  are  many  commercial  filter  design  packages  available  in  the  market,  for  example.  Signal 
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Processing  Tool  Box  in  MATLAB  (PROMATLAB  1991).  Low-pass  or  band-pass  filters  are  often  used 
to  eliminate  noise.  In  general,  the  type  of  filter  used  is  not  critical.  The  important  parameter  is  the  proper 
choice  of  the  break  fr^uency  of  the  filters,  since  if  the  input  is  not  a  band  limited  signal,  the  filter  will 
suppress  useful  information  at  the  same  time  as  the  noise  is  reduced.  The  filter  must  reduce  the  effects 
of  the  external  disturbances  and  those  of  unmodeled  dynamics,  thereby  preventing  sudden  variations  in 
the  estimated  parameters.  As  a  rule  of  thumb,  the  lower  break  frequency  should  be  at  least  1  decade 
below  the  desired  closed  loop  gain  crossover  ftequency  (the  frequency  at  which  the  amplitude  ratio  of  the 
control  loops  is  zero  db),  and  further,  the  upper  break  frequency  should  be  about  2  to  10  times  the  gain 
crossover  frequency.  If  is  too  high,  the  estimator  may  still  attempt  to  fit  the  model  to  high  frequency 
content  noise  and  plant  dynamics. 


Normalization 

In  addition  to  slow  sampling  and  prefiltering  methods  mentioned  above,  there  are  some  other 
methods,  such  as  normalization,  data  scaling,  dead-zone,  projection,  etc.,  to  robustify  the  identifier  when 
noise  and  urunodeled  dynamics  are  present.  It  is  well  known  that,  in  practical  applications  of  least-squares 
estimation,  poor  data  often  lead  to  numerical  problems.  This  is  especially  true  if  the  excitation  of  the 
process  is  poor.  Normalization  methods  are  used  in  this  algorithm  to  improve  the  numerical  quality  of 
the  computations  and  provide  better  estimation  ccmvergence  (Middleton  and  Goodwin  1990;  Sripada  and 
Fisher  1987).  The  normalization  method  normalizes  the  regressor  vector  F(t),  which  is  composed  of  input 
and  output  data  sets  {u(.)}  and  {y(.)},  by  the  normalization  factor  n(t)  =  max(l,  IIF(t)ll)  so  that  all  the 
elemerts  of  F„(t)  =  F(t)/n(t)  <  1.  The  normalization  reduces  the  problem  of  unbounded  disturbances  to  a 
problem  of  bounded  disturbances. 


Time  Delay  and  Model  Order 

Determination  of  the  correct  time  delay  is  an  important  issue  in  the  modeling  procedure.  For  boiler 
systems,  time  delay  usually  varies  with  different  operating  regimes.  A  fixed  time  delay  for  the  whole 
operating  range  is  not  possible.  One  solution  to  this  problem  is  the  use  of  a  higher  order  model  to  allow 
accommodation  to  varying  time  delay.  For  example,  the  water  level  loop  basically  behaves  as  an 
integrator  with  nonminimum  phase  response,  so  a  first-order  model  could  be  used.  To  accommodate  the 
varying  time  delay,  a  third-order  model  was  used  to  model  this  loop.  Simulation  results  show  that  with 
such  a  model  the  system  could  adapt  to  a  fixed  structural  time  delay  varying  from  10  to  23  seconds.  But 
if  the  time  delay  were  overestimated  further,  the  step  response  of  this  loop  will  become  unreasonable,  such 
as  when  die  time  delay  was  set  equal  to  24  seconds.  This  suggests  that  when  choosing  a  time  delay,  it 
is  better  to  choose  too  short  a  delay  rather  than  too  long  a  delay.  A  higher  order  model  must  be  used  to 
increase  the  flexibility  so  that  a  varying  time  delay  can  be  accommodated.  In  our  simulation,  when  a 
lower  order  model  was  used,  it  seems  that  the  plant  still  would  be  able  to  accommodate  the  time  delay 
varying  from  10  to  23  seconds. 
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7  FAULT  TOLERANCE  AND  DIAGNOSTICS 


The  purpose  of  the  fault  tolerance  system  is  to  provide  gracefully  degrading  performance  under  the 
deterioration  and  failure  of  both  hardware  and  software  system  components,  an^or  incorrect  operator 
acticms,  and  under  such  conditions,  to  indicate  the  faulty  component.  Here  we  only  consider  the  simplest 
fault  tolerant  features  that  could  be  implemented  in  DMACS  supervisory  software. 

Typical  failure  modes  for  boiler  control  systems  include: 

1.  Actuator  failure,  e.g.,  actuator  sticking,  motor  burnout,  and  diaphragm  leakage  in  pneumatic 
actuators. 

2.  Sensor  (transmitter)  failure  or  corruption  of  sensor  (pressure,water  level,  oxygen)  reading. 

3.  Deterioration  of  fuel  quality,  e.g.,  dilution  of  natural  gas  with  additives,  dilution  of  oil  with 
water  in  oil-fired  boilers,  burner  clogging. 

4.  Refractory  breakdown  that  exposes  boiler  walls  to  damaging  temperatures. 

5.  Leakage  of  pressurized  water  or  steam  from  sight  glass  or  superheater  tubes. 

6.  Controller  failure  (software  and/or  hardware). 

Procedures  listed  below  in  the  order  corresponding  to  failure  modes  given  above  addressed  each 
failure  type. 

1.  Actuator  deterioration  is  addressed  by  closing  the  local  loop  around  each  actuator  (implemented 
in  pMac  code),  which  ensures  faidiful  execution  of  the  control  signal  under  significant  changes  in  actuator 
characteristics.  The  performance  of  these  local  loops  is  continually  monitored,  and  when  the  performance 
measure  (e.g.,  ITAE  criterion,  the  Integral  of  Time  Multiplied  by  Absolute  Error  Criterion  (Dorf  1989): 

T 

I  =  Jtle(t)|dt  [Eq  56] 

0 


exceeds  a  chosen  threshold,  the  actuator-specific  failure  signal  is  displayed  on  screen  indicating  that 
actuator  performance  is  unacceptable. 

2.  Sensor  failure  is  addressed  by  checking  if  sensor  readings  stay  within  the  physically  meaningful 
bounds.  If  bounds  are  exceeded,  the  sensor-specific  failure  signal  is  displayed.  The  graceful  degradation 
of  sensor  readings  could  be  carried  out  by  combining  the  main  sensor  data  with  the  data  from  a  redundant 
set  of  sensors  to  arrive  at  compromise  sensor  reading.  A  plant  model  could  be  used  to  predict  ouqtuts  for 
knowm  inputs  and  compare  them  to  actual  output  readings  for  discrejiancy.  If  the  discrepancy  exceeds 
a  threshold,  a  failure  signal  is  displayed. 

3.  Fuel  quality  and  burner  clogging  are  monitored  by  periodically  comparing  normal  projected 
boiler  efficiency  with  the  observed  efficiency  both  computed  on-line.  A  second  technique  involves 
checking  for  excessive  fuel  demand  for  the  current  setpoint. 
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4.  Refractory  breakdown  can  be  detected  by  the  discrepancy  between  nominal  and  actual  air/fuel 
demand  for  the  maintenance  of  a  particular  setpoint 

5.  Leakage  is  detected  by  the  pressure  loss  under  nominal  or  excessive  fiiel/air  consumption  for 
a  particular  setpoint. 

6.  Controller  failure  is  determined  by  comparing  ccmtroller  ou^ut  with  baseline  (normal  operaticm) 
for  a  given  input  and  setpoint  using  some  (e.g.,  ITAE)  performance  criterion.  If  an  adaptive  algorithm 
fails,  control  is  switched  to  a  PID  ccMitroller  in  DMACS,  which  provides  acceptable  performance. 

Identifiers  in  the  pMac  code  provide  information  for  pattern  recognition,  which  could  map  parameter 
changes  into  a  variety  of  failure  modes  and  requests  for  maintenance  for  various  subsystems.  This  aspect 
should  be  explored  further. 
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8  CONCLUSIONS  AND  RECOMMENDATIONS 


Hiis  study  developed  a  prototype  Central  Energy  Plant  Adaptive  Control  System  using  a 
microconputer  in  conjunction  with  thiid-party  control  software  to  collect  data,  perfonn  advanced  control 
ftmctions,  generate  reports,  and  to  serve  as  an  alarm  and  an  operator  interface. 

The  prototype  boiler  control  system  was  successfully  field  tested  in  July  1991  at  the  University  of 
Illinois  AU)bot  Power  Plant.  The  system  demonstrated  its  ability  to  control  steam  pressure  and  stack 
oxygen  levels  significantly  better  than  that  of  the  conventional  PID  system.  While  the  water  level  control 
loop  and  the  excess  oxygen  loop  showed  deficiencies  during  abrupt  steam  demand  reduction,  these 
occurrences  may  have  been  caused  by  the  slow  sampling  rate  used  during  the  test  period.  Identification 
protection  schemes  were  implemented  to  permit  higher  sampling  rates,  and  it  is  anticipated  that  the 
reconfigured  controller  will  perform  well  under  all  conditions  (Chapter  4). 

During  the  field  test,  the  controller  hardware  was  critically  examined  and  improvements  were 
identified  and  made.  Automatic  startup  and  shutdown  procedures  were  implemented,  and  possible  self¬ 
diagnostic  techniques  were  explored  (Chapter  3). 

The  great  advantage  of  the  GPC  system  lies  in  its  control  loops,  which  exhibit  significant  non- 
linearities  over  the  entire  operating  range.  It  is  recommended  that  a  longer  term  field  demonstration  be 
done  to  apply  GPC  techniques  to  an  operating  unit  Such  a  dononstration  would  ideally  be  perfonned 
in  a  location  staffed  by  experienced  operators  arul  maintenance  personnel.  The  project  should  exploit  the 
knowledge  and  experience  of  academic  developers,  control  system  manufacturers,  and  a  consulting 
engineer  to  develop  a  practical  demonstration  projea  of  this  promising  technology. 
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APPENDIX  A:  Plan  for  July  Experiment:  Closed  Loop  Control  System  Testing  on  Abbott  Boiler 
No.  2  (Revised  9  July  1991) 


Overall  Experimental  Procedure  (GPC  Test — Wednesday,  10  July) 

1.  Operators  should  set  the  oxygen  for  boiler  No.  2  at  a  moderately  high  value,  and  stabilize  the 
drum  level  with  moderate  boiler  load.  Other  boilers  on  the  header  may  be  on  automatic  during  the  inner 
loop  tuning  imx:edure,  but  all  M/A  stations  for  boiler  No.  2  should  be  on  manual  at  this  time. 

2.  Record  the  current  time  (use  the  AT’s  clock)  and  start  the  boiler  control  program  (GPC  version) 
on  the  Micromac. 

3.  Check  the  PC’s  hard  disks  for  sufficient  space  to  store  data. 

4.  Startup  the  FIX  on  the  PC. 

a.  Start  historical  trend  collect  ruiuiing  in  the  background. 

b.  Enter  view  GPC  1. 

5.  Execute  a  bumpless  transfer  to  micromac  control  (Check  that  all  M/A  stations  are  on  manual. 
Switch  the  selector  to  “Mac.”  Match  the  command  signals  and  switch  the  gas,  air  and  water  flow  M/A 
stations  to  “Auto”). 

6.  Execute  the  Inner  Loop  Startup  Procedure  detailed  below.  At  this  point,  the  actuator  controllers 
are  in  place  on  each  subsystem. 

7.  Execute  the  Outer  Loop  Startup  Procedure  (GPC  case)  detailed  below.  At  this  point,  both  the 
actuatOT  controllers  and  process  controllers  are  in  plia%  cm  each  subsystem. 

8.  Record  the  parameter  values  for  each  loop. 

9.  Turn  on  the  PC  data  log. 

10.  Put  each  subsystem  through  a  sequence  of  small  process  setpoint  changes  to  record  its  perfor¬ 
mance. 


11.  Apply  a  sequence  of  small  steam  demand  disturbances  by  changing  the  fuel  setting  cm  boiler 
No.  3.  Record  the  fuel  settings. 

12.  Return  each  subsystem  to  mode  4. 

13.  Turn  off  the  PC  data  log  and  save  the  file  obtained  Gradually  shift  the  distribution  of  load  by 
changing  die  fuel  setting  cm  boiler  No.  3  so  as  to  bring  boiler  No.  2  into  a  different  operating  regime. 
Manually  cxmtrol  actuator  setpoints  from  Fix  to  maintain  pressure,  oxygen  and  drum  level  near  desired 
values. 

14.  Reset  the  covariances  cm  the  three  cxiter  loops. 

15.  Re-execute  the  Outer  Lcxip  Startup  Pixx;edure  (GPC  case)  detailed  below.  At  this  point,  both  the 
actuator  controllers  and  process  controllers  are  in  place  on  each  subsystem. 
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16.  Record  the  parameter  values  for  each  loop. 

17.  Turn  cm  the  PC  data  log. 

18.  Put  each  subsystem  through  a  sequence  of  small  process  setpoint  changes  to  record  its  tracking 
performance. 

19.  Apply  a  sequence  of  small  steam  demand  disturbances  by  changing  the  fuel  setting  on  boiler 
No.  3.  Reco^  the  fuel  settings. 

20.  Gradually  return  die  steam  load  back  to  its  original  distribution  while  leaving  the  outer  loop 
controllers  at  their  current  settings,  note  whether  good  performance  is  maintained. 

21.  Put  each  subsystem  through  a  sequence  of  small  process  setpoint  changes  to  record  its  tracking 
performance. 

22.  Apply  a  sequence  of  small  steam  demand  disturbances  by  changing  the  fuel  to  record  its 
tracking  performance. 


Overall  Experimental  Procedure  (PID  Test — ^Thursday,  11  July) 

1.  Operators  should  set  the  oxygen  for  boiler  No.  2  at  a  moderately  high  value,  and  stabilize  the 
drum  level  with  load  chi  boiler  No.  2  the  same  as  for  the  GPC  test.  Other  boilers  on  the  header  may  be 
in  automatic  during  the  inner  Icxip  tuning  procedure,  but  all  M/A  staticms  for  boiler  No.  2  shcxild  be  on 
manual  at  this  time. 

2.  Record  the  current  time  (use  the  AT’s  clock)  and  start  the  boiler  control  program  (PID  versicHi) 
cm  the  Micromac. 

3.  CTheck  the  PC’s  hard  disks  for  sufficient  space  to  store  data. 

4.  Startup  die  FIX  cm  the  PC. 

a.  Start  historical  trend  collect  running  in  the  background. 

b.  &iter  view  PIDl. 

5.  Execute  a  bumpless  transfer  to  micromac  ccHitrol  (Check  that  all  M/A  staticms  are  in  manual. 
Switch  the  selector  to  "Mac.”  Match  the  command  signals  and  switch  the  gas,  air  and  water  flow  M/A 
stations  to  “Auto”), 

6.  Execute  the  Inner  Lcmp  Startup  Procedure  detailed  below.  At  this  point,  die  actuator  ccmtrollers 
are  in  place  on  each  subsystem. 

7.  Execute  the  Outer  Loop  Startup  Procedure  (PID  case)  detailed  below.  At  this  point,  both  the 
actuator  controllers  and  process  controllers  are  in  place  on  each  subsystem. 

8.  Turn  cm  die  PC  data  log. 
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9.  Put  each  subsystem  through  the  same  sequence  of  small  fuocess  setpoint  changes  as  was  used 
in  the  GPC  experiment. 

10.  A{^ly  a  sequence  of  small  steam  demand  disturbances  by  changing  the  fuel  setting  on  boiler 
No.  3  in  the  same  sequence  as  was  used  in  the  GPC  experiment  (if  possible). 

1 1.  Return  each  subsystem  to  mode  4. 

12.  Turn  off  the  PC  data  log  aiul  save  die  file  obtained. 

13.  Gradually  diift  the  distribution  of  load  by  changing  the  fuel  setting  on  boiler  No.  3  in  die  sune 
sequence  as  was  used  in  the  GPC  experiment  (if  possible)  so  as  to  bring  boiler  No.  2  into  a  different 
operating  regime  (if  possible,  to  the  same  regime  as  was  used  in  the  GPC  experiment).  Manually  control 
the  actuator  setpoints  from  the  FIX  to  maintain  pressure,  oxygen,  and  drum  level  near  desired  values. 

14.  Re-execute  the  Outer  Loop  Startup  Procedure  (PID  case)  detailed  below.  At  this  point,  both  the 
actuator  ccmtrollers  and  process  controllers  are  in  place  on  e^h  subsystem. 

15.  Turn  on  the  PC  data  log. 

16.  Put  each  subsystem  through  a  sequence  of  small  imicess  setpoint  changes  to  record  its  tracking 
performance. 

17.  Af^ly  a  sequence  of  small  steam  demand  disturbances  by  changing  the  fuel  setting  on  boiler 
No.  3  in  the  same  sequence  as  was  used  in  the  GPC  experimoit  (if  possible). 

18.  Gradually  return  the  steam  load  back  to  its  original  distribution  while  leaving  the  outer  loop 
controllers  at  their  current  settings.  Note  whether  good  performance  is  maintained. 

19.  Put  each  subsystem  through  a  sequence  of  stiudl  process  setpoint  changes  to  record  its  tracking 
performance. 

20.  Apply  a  sequence  of  small  steam  denumd  disturbances  by  changing  the  fuel  setting  on  boiler 
No.  3  in  die  same  sequence  as  was  used  in  die  GPC  experiment  (if  possible). 


Inner  Loop  Startup  Procedure-Using  PID’s  on  AU  Actuator  Loops 

1.  Run  die  batch  file  “sqnesp”  from  the  Compaq.  The  conpaq  will  begin  logging  all  I/O  data  via 
Kermit 

2.  Step  the  first  command  (fuel  flow,  air  flow,  or  water  flow)  by  4  percent  in  a  safe  I/O  via 
Kmnit. 

3.  Wait  for  the  associated  measured  inocess  input  flow  to  settle. 

4.  Stq>  the  command  back  to  its  original  valtw. 

5.  Stop  the  data  logging  on  Compaq.  This  will  automatically  start  Matlab,  which  runs  an  m-file 
and  permits  a  view  of  die  plots  for  delay  arid  maximum  slqie,  and  conqnites  initial  estimates  for  die  PID 
parameters  via  die  Ziegler-Nichols  rules. 
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6.  G>inpare  with  plots  and  gains  obtained  off-line. 

7.  Exit  matlab  on  the  Compaq. 

8.  If  the  step  response  data  is  satisfactory,  get  a  hard  copy: 

a.  “print  sqnesp” 

b.  “graph  matl^”. 

9.  Save  die  step  response  data,  e.g.,  “copy  sq)resp.mat  dataMfstepl.mat”. 

10.  At  the  FIX  keyboard,  enter  initial  PID  gains — either  from  step  response  program  or  off-line 
tuning,  but  cut  the  overall  gain  by  a  factor  of  10  to  ensure  stability. 

1 1.  Turn  on  the  actuator  controller  (switch  to  mode  3 — inner  loop  control  mode). 

12.  When  the  measured  process  input  reaches  steady  state,  run  the  batch  file  “turtresp”  from  the 
Compaq.  The  Compaq  will  begin  logging  all  I/O  data  via  Kermit 

13.  Step  the  first  setpoint  by  S  percent  in  a  safe  directimi. 

14.  Wait  for  the  associated  measured  process  input  flow  to  settle. 

15.  Step  the  setpoint  back  to  its  original  value. 

16.  Stop  the  data  logging  on  the  Compaq.  This  will  automatically  start  Matlab,  which  runs  an  m-file 
and  permits  a  view  of  a  plot  of  the  closed-loop  response,  and  computes  the  rise  time  and  percent 
overshoot 

17.  Compare  with  plots  and  performance  measures  obtained  off-line. 

18.  Exit  matlab  on  the  Compaq. 

19.  If  the  actuator  response  data  is  satisfactory,  get  a  hard  copy: 

a.  “print  actresp” 

b,  “grsqrh  matlab”. 

20.  Save  the  actuator  response  data,  e.g.,  “copy  actresp.mat  dataVffrespl.mat”. 

21.  From  the  FIX  keyboard,  adjust  the  PID  gains. 

22.  When  the  measured  process  input  again  reaches  steady  state,  re-run  the  batch  file  “actresp”  from 
the  Compaq.  The  cmnpaq  will  again  begin  logging  all  I/O  data  via  Kermit.  Repeat  steps  17  to  26  until 
the  clos^-loop  response  is  satisfactory. 

23.  Mtttch  die  current  command  with  die  ctmtroUer  output  and  switch  die  current  loop  back  into 
mode  2. 

24.  Repeat  steps  6  to  28  for  each  inner  loop  (fuel  flow,  air  flow,  water  flow). 
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Older  Loqp  Procedure  (GPC  Case) 


Begin  in  with  all  subsystems  in  mode  3 — inner  loop  control  mode.  At  this  time,  boiler  No.  3  is  still 
in  automatic  (i.e.,  maintaining  pressure). 

1.  Put  boiler  No.  3  on  manual. 

2.  At  the  FIX,  enter  view  GPC2  via  the  PgDn  key. 

3.  Switch  all  subsystems  to  mode  4 — outer  loop  identification  mode. 

4.  Select  the  current  loop  to  tune  (e.g.,  drum  level — loop  S). 

5.  Run  the  step  response  batch  file  for  the  cunent  loop  (e.g.,  “wistep”)  finran  the  Compaq.  The 
Compaq  will  begin  logging  all  I/O  data  via  Kermit. 

6.  Step  die  actuator  setpoint  for  the  current  loop  (e.g.,  water  flow  setpoint)  by  4  peicoit  in  a  safe 
directioiL 

7.  Wait  for  the  associated  measured  process  output  (drum  pressure,  oxygen,  drum  level)  to  settle. 

8.  Step  the  actuator  setpoint  back  to  its  original  value. 

9.  Stop  the  data  logging  on  the  Compaq.  This  will  automatically  start  Matlab,  which  runs  an  m-flle 
that  permits  us  to  view  dK  {dots  for  delay.  Check  the  measured  delay  against  the  delay  assumed  in  the 
boiler  control  program.  Change  the  assumed  delay  if  necessary. 

10.  Compare  with  step  response  plots  obtained  off-line. 

1 1.  Exit  matlab  on  the  Compaq. 

12.  If  the  stq)  response  data  is  satisfactory,  get  a  hard  copy: 

a.  e.g.,  “print  wistep” 

b.  “graph  matlab”. 

13.  Save  the  step  re^nse  data,  e.g.,  “copy  wlstep.mat  dataXwlstepl.mat”. 

14.  Provide  further  excitation  to  actuator  setpoints. 

15.  When  the  parameters  for  a  loop  converge  (i.e.,  parameter  values  appear  stable  and  give 
reastmalde  gains  and  rise  times,  prediction  error  is  small,  covariance  mean  diagon^  is  small),  turn  on  the 
outer  kK^’s  controller  (i.e.,  switch  to  mode  S — control  mode). 

16.  When  the  measured  process  output  reaches  steady  state,  run  the  associated  closed-loop  response 
batch  file  (e.g.,  “wlresp”)  from  the  Compaq.  The  Compaq  will  begin  logging  all  I/O  data  via  ^rmit 

17.  Step  the  process  setpoint  by  a  reasonable  amount  (e.g.,  5  psi  for  drum  pressure,  2  percent  for 
oxygen,  OJ  in.  for  drum  level)  in  a  safe  direction. 

18.  Wait  for  the  measured  process  oufout  to  %ttle. 
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19.  Step  the  process  setpoint  back  to  its  oii^nal  value. 

20.  Stop  the  data  logging  on  the  Compaq.  This  will  automatically  start  Matlab,  which  runs  an  m-file 
that  permits  us  to  view  a  plot  of  the  closed-loop  response,  and  computes  the  rise  time  and  percent 
overshoot 

21.  Compare  with  plots  and  performance  measures  obtained  off-lirK. 

22.  Exit  matlab  on  the  Compaq. 

23.  If  the  outer  loop  response  data  is  satisfactory,  get  a  hard  copy; 

a.  e.g.,  “print  wlresp” 

b.  “graph  matlab". 

24.  Save  the  actuator  response  data,  e.g.,  “copy  wlresp.mat  dataXwlrespl.mat". 

25.  If  the  response  is  unsatisfactory,  return  to  mode  4  and  repeat  steps  14  to  23. 

26.  If  the  response  is  satisfactory,  leave  die  current  loop  in  mode  5. 

27.  Repeat  the  above  procedure  for  the  next  outer  loop  (drum  water  level,  drum  pressure,  oxygen). 
Outer  Loop  Startup  Procedure  (PID  Case) 

Begin  with  all  subsystems  in  mode  3 — inner  loop  control  mode.  At  this  time,  boiler  No.  3  is  still 
on  automatic  (i.e.,  maintaining  pressure). 

1.  Put  boiler  No.  3  on  manual. 

2.  At  the  FIX,  enter  view  PID2  via  die  PgDn  key. 

3.  Switch  all  subsystems  to  mode  4 — outer  loop  identiflcatiCHi  mode. 

4.  Select  die  current  loop  to  tune  (e.g.,  drum  level — ^loop  S). 

5.  Run  die  step  response  batch  file  for  die  current  loop  (e.g.,  “wlstep")  from  the  Compaq.  The 
Compaq  will  begin  logging  all  I/O  data  via  Kermit 

6.  Step  first  actuator  setpoint  (fuel  flow,  air  flow,  or  water  flow)  by  4  percent  in  a  safe  direction. 

7.  Wait  for  the  associated  measured  process  output  (drum  pressure,  oxygen,  drum  level)  to  settle. 

8.  Step  the  actuator  setpoint  back  to  its  original  value. 

9.  Stop  the  data  logging  on  Compaq.  This  wall  automatically  start  Matlab,  which  runs  an  m-file 
diat  permits  us  to  view  the  plots  for  delay  a^  maximum  slope,  and  computes  initial  estimates  for  the  PID 
param^rs  via  the  Ziegler-Nichols  rules. 
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10.  Compare  with  plots  and  gains  obtained  off-line. 

11.  Exit  matlab  on  the  Compaq. 

12.  If  the  step  response  data  is  satisfactory,  get  a  hard  copy: 

a.  e.g.,  “print  wlstep” 

b.  “graph  matlab". 

13.  Save  the  step  response  data,  e.g.,  “copy  wlstep.mat  dataNwlstepl.maL” 

14.  At  the  FIX  keyboard,  enter  initial  PID  gains — either  from  step  re^nse  program  or  off-line 
tuning,  but  cut  the  overall  gain  by  a  factor  of  10  to  ensure  stability. 

15.  Turn  on  the  outer-loop  controller  (switch  to  mode  5 — full  control  mode). 

16.  When  the  measured  process  output  reaches  steady  state,  run  the  associated  closed-loop  response 
batch  file  (e.g.,  “wlresp”)  from  the  Compaq.  The  Compaq  will  begin  logging  all  lAD  data  via  KenniL 

17.  Step  the  first  setpoint  by  a  reasonable  amount  (e.g.,  5  psi  for  drum  pressure.  2  percent  for 
oxygen,  OJ  ia  for  drum  level)  in  a  safe  directioa 

18.  Wait  for  the  measured  process  output  to  settle. 

19.  Step  the  setpoint  back  to  its  original  value. 

20.  Stop  the  data  logging  on  the  Compaq.  This  will  automatically  start  Matlab,  which  runs  an  m-file 
that  permits  us  to  view  a  plot  of  the  closed-loop  response,  and  computes  the  rise  time  and  percent 
overdnot 

21.  (i;ompare  with  plots  and  performance  measures  obtained  off-line. 

22.  Exit  matlab  on  the  Compaq. 

23.  If  the  outer  loop  response  data  is  satisfactory,  get  a  hard  copy: 

a.  e.g.,  “print  wlresp" 

b.  “graph  matlab”. 

24.  Save  the  actuator  response  data,  e.g.,  “copy  wlresp.mat  dataXwlrespl.mat”. 

25.  From  the  FIX  keyboard,  adjust  the  PID  gaits. 

26.  When  the  measured  process  itqnit  again  reaches  steady  state,  re-r\m  the  batch  file  “actresp”  from 
the  cmnpaq.  The  Compaq  will  again  begin  logging  the  loop’s  I/O  data  via  Kermit  Repeat  steps  17  to  26 
urttil  die  closed-loop  respmse  is  satisfactory. 

27.  Leave  the  current  loop  back  in  mode  5. 

28.  Rqieat  steps  6  to  28  for  the  next  outer  loop  (drum  pressure,  oxygen,  drum  level). 
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APPENDIX  B:  PID-ControUcr  Parameters  Tuning  Rules 

In  general,  there  are  two  basic  methods  for  tuning  PID-controller  parameters:  the  transient-respcmse 
method  uid  the  ultimative-sensitivity  method  (Astrbm  and  Wittenmarii  1984).  The  continucHis-time  PID 
cmitroller  is  often  written  in  Laplace  form  as; 


U(s)  =  K(l 


V 

1  +TpS./N 


)E(s) 


[EqBl] 


Where: 

lUs)  and  E(s)  are  the  Laplace  transforms  of  die  controller  output  and  the  error  signal,  respectively 
K  is  the  proportional  gain 
T,  is  the  integral  time 
Tp  is  the  derivative  time. 

In  the  controller  there  is  a  filter,  with  time  constant  for  the  derivative  part.  Also,  constant  N 
is  often  in  the  range  3  to  10. 

The  digital  PID-controller  is  parameterized  by  approximating  the  derivative  with  a  forward 
difference  (Euler’s  method)  as: 

u(KT)  =  KJl  ^  T  •.  ,^“^t:,!!.]e(KT)  [Eq  B2] 

Ty(q-1)  T(q+Y) 


Where: 

u(kT)  and  e(kT)  are  the  controller  output  and  the  error  signal  at  time  kT 
it  is  the  sampling  numbers 
T  is  the  sampling  period 
Kj  is  the  proportional  gain 

Tjd  and  are  the  discrete-time  equivalent  to  integral  and  derivative  times 
9  is  a  difference  operator 
■p-expf-TNA’d)  is  a  constant. 


The  Transient-Response  Method 

To  use  the  transient-response  method,  the  steepest  slope,  R,  and  the  delay  time,  L,  are  measured 
from  a  unit-step  response  of  the  open-loop  system  (Figure  B 1 ).  The  parameters  for  the  PID-controller  are  ^ 
then  obtained  from  Table  Bl. 
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Figure  Bl.  Measurement  on  the  Open  Loop  Untt*Step  Response  of  the  Variables  R  and  L  for  the 
Transient-Response  Method. 


Table  Bl 

Controller  Parameters  When  Using 
the  Transient-Response  Method 


Kd 

Tid 

Tdd 

p 

1/RL 

PI 

0.9/RL 

3L 

PID 

1.2/RL 

2L 

0.5L 

The  Ultimatiye-Sensitiyity  Method 

In  the  ultimative-sensitivity  method,  a  P-c<HitroiIer  is  used  first  to  control  die  system.  The  gain  of 
the  controller,  and  the  period  time,  Tp,  are  measured  when  the  closed-loop  system  is  on  the  stability 
boundary.  The  parameters  of  the  controller  are  then  obtained  from  Table  B2. 


Table  B2 

Controller  Parameters  When  Using  the  UHimative-Sensitivity  Method 


Kd 

Tid 

Tdd 

P 

0JK«. 

PI 

0.45K^ 

y\2 

HD 

0.6K«. 

T^ 

The  tuning  rules  above  should  be  used  only  as  a  first  approximation.  The  final  tuning  usually  has 
to  be  done  manually.  In  general,  the  transient-response  method  is  easy  to  use,  since  it  only  needs  to 
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measure  the  uitit-step  response  of  tfie  open  loop  system.  The  ultimative-sensitivity  method,  however, 
requires  experimentation  on  the  closed  loop  system  to  find  the  maximum  gain  kmax  and  the  oscillatory 
period  tima  Tp .  This  is  often  undesirable  in  the  industrial  processes  for  safety  reasons.  Hence,  the 
ultimative-sensitivity  method  is  only  used  for  processes  with  slow  dynaimcs  such  as  those  that  occur  in 
industrial  furnaces. 
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APPENDIX  C:  Database  fM*  Startup  and  Shutdown  Procedures 


TAG .  SS-MODE 

DESC .  SELECTION  OF  STARTUP  OR  SHUTDOWN 

DV-DEVICE .  AD6 

HT-H/W  OPTIONS..  MAC6000 

ICLADDR .  1:0:145 

SC-SIG  COND.....  NONE 

EL-LO  EGU .  0.00 

EH-HI  EGU .  100.00 

ET-BGU  TAG . 

PA-PLANT  AREA...  ALL 


AE-ALM  ENABLE...  DISABLE 

NX-NEXT  BLK . 

LL-LO  OP  LIMIT..  0.00 
LH-HI  OP  LIMIT..  100.00 
LR-RATE  LIMIT...  100.00 
CS-COLD  START...  0.00 
WS-WARM  START...  NO 
OR-OUT  REVERSE..  NO 


TAG .  PROMPT 

DESC .  STATE  OF  PROMPT 

DV-DEVICE AD6 

HT-H/W  OPTIONS..  MAC6000 

lO-ADDR .  1:0:146 

SC-SlG  COND .  NONE 

EL-LO  EGU .  0.00 

EH-HI  EGU .  100.00 

ET-EGU  TAG . 

PA-PLANT  AREA...  ALL 
AE-ALM  ENABLE..  DISABLE 

NX-NEXT  BLK . 

LL-LO  OP  LIMIT..  0.00 
LH-HI  OP  LIMIT..  100.00 
LR-RATE  LIMIT...  100.00 
CS-COLD  START...  0.00 
WS-WARM  START...  NO 
OR-OUT  REVERSE.  NO 


TAG .  FLAG 

DESC .  PROMPT  FLAG 

DV-DEVICE. .  AD6 

HT-H/W  OPTIONS..  MAC6000 

lO-ADDR .  1:0:147 

SC-SKj  COND.....  NONE 
EL-LO  EGU .  0.00 


miT 


mrr 


ET-EGU  TAG . 

PA-PLANT  AREA...  ALL 
AE-ALM  ENABLE...  DISABLE 

NX-NEXT  BLK . 

LL-LO  OP  LIMIT..  0.00 
LH-m  OP  LIMIT..  100.00 
LR-RATE  LIMIT...  100.00 
CS-COLD  START...  5.00 
WS-WARM  START...  NO 
OR-OUT  REVERSE..  NO 


TAG .  A.PROMPT 

DESC .  ACHO  VALUE  OF  PROMPT  FROM  UMAC 

DV-DEVICE .  AD6 

HT-H/W  OPTIONS..  MAC6000 

lO-ADDR .  1:0:148 

SC-SIG  CX)ND .  NONE 

EL-LO  ECU .  0.00 

EH-HI  EGU .  100.00 


ET-EGU  TAG . 

SM-SMOOTHING....  0 
IS-INTT  SCAN....  ON 
ST-SCAN  TIME....  5 

NX-NEXT  BLK . 

LL-LO  LO  ALM....  0.00 

AL-LO  ALARM .  0.00 

AH-m  ALARM.....  100.00 
HH-m  HI  ALM....  100.00 


RC-ROC  ALM .  0.00 

DB-1£AD  BAND....  0.20 
PA-PLANT  AREA...  ALL 

AP-ALM  PRI .  L 

AE-ALM  ENABLE...  ENABLE 
lA-INTT  A/M.....  AUTO 

TAG .  A_FLAG 

DESC .  Acho  value  of  FLAG  from  umac 

DV-DEVICE .  AD6 

HTT-H/W  OPTIONS..  MAC6000 

lO-ADDR. .  1:0:149 

SC-SIG  COND.....  NONE 

EL-LO  EGU .  0.00 

EH-m  EGU.......  100.00 

ET-EGU  TAG . 


SM-SMOOTHING....  0 
is-wrr  SCAN....  on 
ST-SCAN  TIME....  5 
NX-NEXT  BLK . 
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LL>LO  LO  ALM.... 

0.00 

AL-LO  ALARM . 

0.00 

AH-m  ALARM..... 

100.00 

HH-Ifl  HI  ALM.... 

100.00 

RC-ROC  ALM . 

0.00 

DB-DEAD  BAND... 

.  0.20 

PA-PLANT  AREA.., 

,  ALL 

AP-ALM  PRI . 

L 

AE-ALM  ENABLE. 

ENABLE 

lA-INTT  A/M . 

AUTO 

TAG . 

MAIN 

PA-PLANT  AREA. 

ALL 

AE-ALM  ENAB... 

ENABLE 

AP-ALM  PRI.... 

L 

ST-SCAN  TIME.. 

5 

B-INIT  SCAN.. 

ON 

lA-INTT  A/M... 

AUTO 

DESC. 

MAIN  PROCEDURES  FOR  STARTUP  &  SHUTDOWN 

00. 

SETUM  0.2000 

01. 

IF  SS-MODE  =  1.00  GOTO  4 

02. 

IF  SS-MODE  =  2.00  GOTO  10 

03. 

GOTO  1 

04. 

CALL  STARTUP 

05. 

SETOUT  PROMPT  28.00 

06. 

DELAY  10 

07. 

SETOUT  SS-MODE  0.00 

08. 

GOTO  1 

09. 

NUL 

10. 

CALLSHUTDN 

11. 

END 

12. 

NUL 

13. 

NUL 

14. 

NUL 

15. 

NUL 

16. 

NUL 

17. 

NUL 

18. 

NUL 

19. 

NUL 

TAG . 

STARTUP 

PA-PLANT  AREA. 

ALL 

AE-ALM  ENAB... 

ENABLE 

AP-ALM  PRI.... 

L 

ST-SCAN  TIME.. 

1 

IS-INTT  SCAN.. 

OFF 

lA-INTT  A/M... 

AUTO 

DESC.. 

STARTUP  PROCEDURES 

00. 

CALL  STARTl 
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01. 

SETOUT  PROMPT  5.00 

02. 

DELAY  10 

03. 

WATTFOR  A_PROMPT  =  0.00 

04. 

CALLSTART2 

05. 

CALLSTART3 

06. 

CALL  START4 

07. 

CALL  START5 

08. 

CALL  START6 

09. 

NUL 

10. 

NUL 

11. 

NUL 

12. 

NUL 

13. 

NUL 

14. 

NUL 

15. 

NUL 

16. 

NUL 

17. 

NUL 

18. 

NUL 

19. 

NUL 

TAG _ 

STARTl 

PA-PLANT  AREA. 

ALL 

AE-ALM  ENAB... 

ENABLE 

AP-ALM  PRI.... 

L 

ST-SCAN  TIME.. 

1 

IS-INTT  SCAN.. 

OFF 

lA-INTT  A/M... 

AUTO 

DESC.. 

SETUP  WF.  WL  &  AF  ACTUATORS 

00. 

SETOUT  WF-SET  50.00 

01. 

SETOUT  WL-SET  0.00 

02. 

SETUM  10.0000 

03. 

WATTFOR  WF  >=  50.00 

04. 

SETLIM  0.2000 

05. 

SETOUT  PROMPT  1.00 

06. 

DELAY  10 

07. 

WATTFOR  A_PROMPT  =  0.00 

08. 

SETOUT  PROMPT  2.00 

09. 

DELAY  10 

10. 

WATTFOR  A.PROMFT  =  0.00 

11. 

SETOUT  AF-SET  40.00 

12. 

SETOUT  PROMPT  3.00 

13. 

DELAY  10 

14. 

WATTFOR  A_PROMPT  =  0.00 

15. 

SETOUT  PROMPT  4.00 

16. 

DELAY  10 

17. 

WATTFOR  A_PROMPT  =  0.00 

18. 

NUL 

19. 

NUL 
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TAG . 

START2 

PA-PLANT  AREA. 

ALL 

AE-ALM  ENAB... 

ENABLE 

AP-ALM  PRI.... 

L 

ST-SCAN  TIME.. 

1 

IS-INTT  SCAN.. 

OFF 

lA-INTT  A/M... 

AUTO 

DESC.. 

PREPROCEDURE  FOR  IGNITE 

00. 

SETOUT  PROMPT  6.00 

01. 

DELAY  10 

02. 

WATTFOR  A_PROMPT  =  0.00 

03. 

SETOUT  PROMPT  7.00 

04. 

DELAY  10 

05. 

WATTFOR  A_PROMPT  =  0.00 

06. 

SETOUT  PROMPT  8.00 

07. 

DELAY  10 

08. 

WATTFOR  A_PROMPT  =  0.00 

09. 

SETOUT  PROMPT  9.00 

10. 

DELAY  10 

11. 

WATTFOR  A_PROMPT  =  0.00 

12. 

SETOUT  PROMPT  10.00 

13. 

DELAY  10 

14. 

WATTFOR  A_PROMPT  =  0.00 

15. 

SETOUT  PROMPT  11.00 

16. 

DELAY  10 

17. 

WATTFOR  A_PROMPT  s  0.00 

18. 

NUL 

19. 

NUL 

TAG . 

START3 

PA-PLANT  AREA. 

ALL 

AE-ALM  ENAB... 

ENABLE 

AP-ALM  PRI.... 

L 

ST-SCAN  TIME.. 

1 

IS-INTT  SCAN.. 

OFF 

lA-INTT  A/M... 

AUTO 

DESC.. 

ATTEMPT  TO  IGNITE  THE  IGNTTOR 

00. 

SETOUT  PROMPT  12.00 

01. 

DELAY  10 

02. 

WATTFOR  A_PROMPT  =  0.00 

03. 

DELAY  10 

04. 

IF  A_FLAG  =  1.00  GOTO  18 

05. 

DELAY  60 

06. 

SETOUT  PROMPT  13.00 

07. 

DELAY  10 

08. 

WATTFOR  A_PROMPT  =  0.00 

09. 

SETOUT  PROMPT  12.00 

10. 

DELAY  10 

11. 

WATTFOR  A_PROMPT  =  0.00 

12. 

DELAY  10 

100 


13. 

14. 

15. 

16. 

17. 

18. 
19. 


IF  A_FLAG  *  1.00  GOTO  18 
SETOUT  PROMPT  14.00 
DELAY  10 

WATTFOR  A_PROMPT  =  0.00 

END 

NUL 

NUL 


TAG . 

PA-PLANT  AREA. 
AE-ALM  ENAB... 
AP-ALM  PRI.... 
ST-SCAN  TIME.. 
IS-INTT  SCAN.. 
lA-INTT  A/M... 
DESC.. 

00. 

01. 

02. 

03. 

04. 

05. 

06. 

07. 

08. 

09. 

10. 

11. 

12. 

13. 

14. 

15. 

16. 

17. 

18. 

19. 


START4 

ALL 

ENABLE 

L 

1 

OFF 

AUTO 

ENGAGE  FF  CONTROL 
SETOUT  PROMPT  15.00 
DELAY  10 

WATTFOR  A_PROMPT  =  0.00 
DELAY  10 

IF  A_FLAG  =  1.00  GOTO  17 
CALL  START2 
CALLSTART3 
SETOUT  PROMPT  15.00 
DELAY  10 

WATTFOR  A.PROMPT  =  0.00 
DELAY  10 

IF  A_FLAG  =  1.00  GOTO  17 
DELAY  60 

SETOUT  PROMPT  16.00 
DELAY  10 

WATTFOR  A_PROMPT  =  0.00 

END 

NUL 

NUL 

NUL 


TAG . 

PA-PLANT  AREA. 
AE-ALM  ENAB... 
AP-ALM  PRI.... 
ST-SCAN  TIME.. 
IS-INTT  SCAN.. 
lA-INTT  A/M... 
DESC.. 

00. 

01. 

02. 

03. 


START5 

ALL 

ENABLE 

L 

1 

OFF 

AUTO 

ENGAGE  FF  CONTROL 
SETOUT  PROMPT  17.00 
DELAY  10 

WATTFOR  A_PROMPT  s  0.00 
SETOUT  PROMPT  18.00 
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04.  DELAY  10 

05.  WATTFOR  A_PROMPT  =  0.00 

06.  SETOUT  FF-SET  FF 

07.  SETOUT  PROMPT  19.00 

08.  DELAY  10 

09.  WATTFOR  A_PROMPT  =  0.00 

10.  SETOUT  PROMPT  20.00 

11.  DELAY  10 

12.  WATTFOR  A_PROMPT  =  0.00 

13.  SETOUT  PROMPT  21.00 

14.  DELAY  10 

15.  WATTFOR  A.PROMPT  =  0.00 

16.  SETOUT  PROMPT  22.00 

17.  DELAY  10 

18.  WATTFOR  A_PROMPT  =  0.00 

19.  NUL 

TAG .  START6 

PA-PLANT  AREA.  ALL 
AE-ALM  ENAB...  ENABLE 
AP-ALM  PRI....  L 
ST-SCAN  TIME..  1 
IS-INTT  SCAN..  OFF 
U-INTT  A/M...  AUTO 

DESC..  ENGADE  OUTER  LOOP  CONTROL 

00.  SETOUT  PROMPT  23.00 

01.  DELAY  10 

02.  WATTFOR  A^PROMPT  =  0.00 

03.  SETOUT  PROMPT  24.00 

04.  DELAY  10 

05.  WATTFOR  A.PROMPT  =  0.00 

06.  SETOUT  M-02-SET  02 

07.  SETOUT  PROMPT  25.00 

08.  DELAY  10 

09.  WATTFOR  A_PROMPT  =  0.00 

10.  SETOUT  OVR-ENGA  0.00 

1 1.  SETOUT  PROMPT  26.00 

12.  DELAY  10 

13.  WATTFOR  A_PROMPT  =  0.00 

14.  SETOUT  DP-SET  DP 

15.  SETOUT  PROMPT  27.00 

16.  DELAY  10 

17.  WATTFOR  AJ*ROMPT  =  0.00 

18.  SETOUT  DP-SET  320.00 

19.  NUL 

TAG .  SHUTDN 

PA-PLANT  AREA.  ALL 
AE-ALM  ENAB...  ENABLE 
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AP-ALM  PRI.... 

L 

ST-SCAN  TIME.. 

1 

IS-INTT  SCAN.. 

OFF 

lA-INTT  A/M... 

AUTO 

DESC.. 

SHUT  DOWN  PROCEDURE 

00. 

NUL 

01. 

CALLSHUTl 

02. 

CALLSHUT2 

03. 

CALL  SHUTS 

04. 

END 

05. 

NUL 

06. 

NUL 

07. 

NUL 

08. 

NUL 

09. 

NUL 

10. 

NUL 

11. 

NUL 

12. 

NUL 

13. 

NUL 

14. 

NUL 

15. 

NUL 

16. 

NUL 

17. 

NUL 

18. 

NUL 

19. 

NUL 

TAG . 

SHUTl 

PA-PLANT  AREA. 

ALL 

AE-ALM  ENAB... 

ENABLE 

AP-ALM  PRI.... 

L 

ST-SCAN  TIME.. 

1 

IS-INTT  SCAN.. 

OFF 

lA-INTT  A/M... 

AUTO 

DESC. 

SHUT  DOWN  DP.  02  &  FF  LOOPS 

00. 

SETOUT  PROMPT  30.00 

01. 

DELAY  10 

02. 

WATTFOR  A_PROMPT  =  0.00 

03. 

SETOUT  DP-SET  0.00 

04. 

DELAY  60 

05. 

SETOUT  PROMPT  31.00 

06. 

DELAY  10 

07. 

WATTFOR  A.PROMPT  =  0.00 

08. 

SETOUT  PROMPT  32.00 

09. 

DELAY  10 

10. 

WATTFOR  A.PROMPT  =  0.00 

11. 

SETOUT  PROMPT  33.00 

12. 

DELAY  10 

13. 

WATTFOR  A_PROMPT  =  0.00 

14. 

SETOUT  FF-SET  0.00 

15. 

DELAY  30 
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16. 

NUL 

17. 

NUL 

18. 

NUL 

19. 

NUL 

TAG . 

SHUT2 

PA-PLANT  AREA. 

ALL 

AE-ALM  ENAB... 

ENABLE 

AP-ALM  PRI.... 

L 

ST-SCAN  TIME.. 

2 

IS-INTT  SCAN.. 

OFF 

lA-INTT  A/M... 

AUTO 

DESC.. 

SHUT  DOWN  THE  REST 

00. 

SETOUT  PROMPT  34.00 

01. 

DELAY  10 

02. 

WATTFOR  A_PROMPT  =  0.00 

03. 

SETOUT  PROMPT  35.00 

04. 

DELAY  10 

05. 

WATTFOR  A_PROMPT  =  0.00 

06. 

SETOUT  PROMPT  36.00 

07. 

DELAY  10 

08. 

WATTFOR  A.PROMPT  =  0.00 

09. 

SETOUT  PROMPT  37.00 

10, 

DELAY  10 

11. 

WATTFOR  A_PROMPT  =  0.00 

12. 

SETOUT  PROMPT  38.00 

13. 

DELAY  10 

14. 

WATTFOR  A_PROMPT  =  0.00 

15. 

NUL 

16. 

NUL 

17. 

NUL 

18. 

NUL 

19. 

NUL 

TAG . 

SHUT3 

PA-PLANT  AREA. 

ALL 

AE-ALM  ENAB... 

ENABLE 

AP-ALM  PRI.... 

L 

ST-SCAN  TIME.. 

1 

IS-INTT  SCAN.. 

OFF 

lA-INTT  A/M... 

AUTO 

DESC.. 

SHUTDOWN  WATER  FLOW 

00. 

SETOUT  PROMPT  39.00 

01. 

DELAY  10 

02. 

WATTFOR  A.PROMPT  =  0.00 

03. 

SETOUT  PROMPT  40.00 

04. 

DELAY  10 

05. 

WATTFOR  A  J»ROMPT  =  0.00 

06. 

SETOUT  WF-SET  0.00 
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07. 

DEI>AY  30 

08. 

SETOUT  PROMPT  41.00 

09. 

DELAY  10 

10. 

WATTFOR  A.PROMPT  =  0.00 

11. 

SETOUT  PROMPT  42.00 

12. 

NUL 

13. 

NUL 

14. 

NUL 

15. 

NUL 

16. 

NUL 

17. 

NUL 

18. 

NUL 

19. 

NUL 
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APPENDIX  O:  FIX  DMACS  Graphic  Environment  for  the  Aut<miatic  Startup  and  Shutdown 
Procedures 


OISP:  PIDl 


ENABLE 


F^re  D2.  Inner  Loop  PID  Control. 
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Figure  D3.  Outer  Loop  PID  Control. 


ENABLE 


OUTPUTS  I  SETPOINTS 


mumm 


DISP:  P1D3 


ENABLE 


Figure  D4.  Output 
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APPENDIX  E:  Prompt  Lists 


PROMPTXrBAD  PROMPT! 

PROMPTOJOB  DONE! 

PROMPTlrPRESS  PgDn  TO  ENGAGE  WF  &  WL  CONTROLLERS 
PROMPT2:OPEN  DAMPER.  OPEN  FLOW  PATH  FROM  AIR  TO  STACK 
PROMPT3:PRESS  PgDn  TO  ENGAGE  AF  CONTROLLER 

PROMPT4:PRECHARGE  THE  BURNER  HEADER  WITH  GAS,  AND  CHECK  FOR  LEAKS 
PROMPT5:VERIFY  THAT  DAMPER  AND  BURNER  REGISTER  ARE  IN  PURGE  POSITION 
PROMPT6:PERFORM  A  UNIT  PURGE 

PROMPT7:CLOSE  THE  MAIN  FUEL  CONTROL  VALVE.  OPEN  THE  SAFETY  SHUTOFF  VALVE 
PROMPT8:ESTABLISH  HEADER  PRESSURE  FOR  BURNER  LIGHTOFF 
PROMPT9:ESTABLISH  IGNITOR  PRESSURE  FOR  IGNITION 
PROMPTlOrSET  AIR  REGISTER  TO  UGHTOFF  POSITION 
PROMPTl  LIGNITE  THE  IGNITOR 

PROMPT12:IF  A  FLAME  IS  NOT  ESTABLISHED.  CLOSE  THE  IGNITOR  SAFETY  VALVE.  ELSE, 
SET  FLAG  =  1 

PROMPT13;RETRY  IGNITE  THE  IGNITOR 

PROMPT14:SE(X>ND  IGNITION  FAILED!  THE  STARTUP  PROCEDURE  WILL  HALT! 
PROMPT15:IF  THE  IGNITOR  IS  OPERATING  CORRECHLY,  SET  FLAG=1  &  SUPPLY  FUEL  TO 
THE  BURNER!  ELSE,  SET  FLAG=0 

PROMPT16:SECOND  FAILURE  TO  UGHTOFF,  THE  STARTUP  PRCX:EDURE  WILL  HALT! 
PROMPT17:HAS  A  STABLE  FLAME  BEEN  ESTABLISHED? 

PROMPT18:SET  AIR  REGISTER  TO  NORMAL  OPERATING  POSITION 
PROMPT19:PRESS  PgDn  TO  ENGAGE  FF  CONTROLLER 
PROMPT20:MAKE  SURE  THAT  IGNITION  IS  NOT  LOST 
PROMPT21:SHUT  DOWN  IGNITOR,  VERIFY  THAT  FLAME  IS  MAINTAINED 
PROMPT22:SHUT  OFF  THE  BURNER  HEADER  ATMOSPHERIC  VALVE 
PROMPT23:HAS  A  STABLE  FLAME  BEEN  ESTABLISHED? 

PROMPT24:RAMP  DOWN  AF  TO  ABOUT  6%  EXCESS  02 
PROMFT2S;PRESS  PgDn  TO  ENGAGE  02  CONTROL  L(X)P 
PROMPT26:IS  THE  BURNER  OPERATING  PROPERLY? 

PROMPT27:PRESS  PgDn  TO  ENGAGE  DP  CONTROL  L(X)P 
PROMPT28:STARTUP  PROCEDURE  DONE! 

PROMPT30:RAMP  DOWN  DRUM  PRESSURE  TO  BRING  THE  BOILER  OFF  LINE. 
PR0MPT31:PRESS  PgDn  TO  DISENGAGE  DP  LOOP  CONTROLLER 
PROMPT32:RAMP  DOWN  FF  Setpoint  UNTIL  AF  REACHES  ITS  PURGE  RATE 
PROMPT33:PRESS  PgDn  TO  DISENGAGE  02  LOOP  CONTROLLER 
PROMPT35:RAMP  DOWN  FF  TO  A  STARTUP  LEVEL 
PROMPT36:PRESS  PgDn  TO  DISENGAGE  FF  CONTROLLER 
PROMPT37:SET  DAMPER  TO  A  STARTUP  POSITION 

PROMPT38:  CLOSE  THE  SAFETY  SHUTOFF  VALVE  TO  SHUT  DOWN  BURNER 
PROMPT39:OPEN  ALL  ATMOSPHERIC  VENT  VALVES 
PROMPT40:COMPLETE  A  UNIT  PURGE 

PR0MPT41:PRESS  PgDn  TO  DISENGAGE  WF  &  WL  CONTROLLERS 
PROMPT42:SHUT  DOWN  BLOWERS  AND  CLOSE  DAMPERS 
PROMPT43:SHUTDOWN  PROCEDURE  DONE! 

SMODEO:WArnNG  INPUT 
SMODELSTARTUP  NOW  ! 

SMOI^:SHUTDOWN  NOW! 
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APPENDIX  F:  Database  for  Effldency  Calculatioo 


TYPE:  AO 

TAG . FUEL 

DV-DEVICE . DUMMY 

ITT-H/W  OPTIONS. 

lO-ADDR . 

SC-SIG  COND.... 

EL-LOEGU .  0 

EH-HIEGU .  5 

ET-EGU  TAG . NO 

OR-OUT  REVERSE.  N 
CS-COLD  START..  2 
DESC....  Fuel  Selection 
LL-LO  OP  LIMIT.  0 
LH-m  OP  LIMIT.  5 
LR-RATE  LIMIT..  5 
NX-NEXT  BLK.... 
CP-COPY  TO . 


TAG . TEMP 

DV-DEVICE. . DUMMY 

HT-H/W  OPTIONS. 

lO-ADDR. . 

SC-SIG  COND.... 

EL-LOEGU .  100.0 

EH-HIEGU .  1000.0 

ET-EGU  TAG . No 

OR-OUT  REVERSE.  N 
CS-COLD  START..  300.0 
IffiSC....  Flue  Gas  Tenqwiatuie 
LL-LO  OP  LIMIT.  100.0 
LH-ffl  OP  LIMIT.  1000.0 
LR-RATE  LIMIT..  900.0 
NX-NEXT  BLK.... 

CP-COPY  TO . 


TAG . 02JUNK 

DV-DEVICE. . DUMMY 

HT-H/W  OPTIONS. 

lO-ADDR. . 

SC-SIG  COND.... 

EL-LOEGU .  0.00 

EH-HIEGU .  20.00 

ET-EGU  TAG . PCT 

QR-OUT  REVERSE.  N 
CS-COLD  START..  5.00 
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DESC....  o2  fm  test 
LL-LO  OP  LIMIT.  0.00 
LH-ffl  OP  LIMIT.  20.00 
LR-RATE  LIMIT..  20.00 
NX-NEXT  BLK.... 
CP-COPY  TO . 


TAG . EFF 

DV-DEVICE. . DUMMY 

HT-H/W  OPTIONS. 

lO-ADDR . 

SC-SIG  CX)ND.... 

EL-LOEGU 0.00 

EH-HIEGU 100.00 

ET-EGU  TAG . PCT 

OR-OUT  REVERSE.  N 
CS-COLD  START.. 

DESC....  Calculated  Combustion  Efficiency 
LL-LO  OP  LIMIT.  0.00 
LH-m  OP  LIMIT.  100.00 
LR-RATE  LIMIT..  100.00 
NX-NEXT  BLK....  EFF-TR 
CP-COPY  TO . 


TAG . C02 

DV-DEVICE . DUMMY 

HT-H/W  OPTIONS. 

lO-ADDR . 

SC-SIG  COND.... 

EL-LOEGU .  0.00 

EH-HIEGU .  100.00 

ET-EGU  TAG . PCT 

OR-OUT  REVERSE.  N 
CS-COLD  START.. 
DESC....  calculated  C02 
LL-LO  OP  LIMIT.  0.00 
LH-HI  OP  LIMIT.  100.00 
LR-RATE  LIMIT..  100.00 
NX-NEXT  BLK.... 
CP-COPY  TO . 


TAG . NG-02 

DV-DEVICE. . DUMMY 

HT-H/W  OPTIONS. 

lO-ADDR. . 

SC-SIG  COND.... 

EL-LOEGU. .  0.00 

EH-HIEGU......  20.00 


no 


ET-EGU  TAG . PCT 

OR-OUT  REVERSE.  N 
CS-COLD  START.. 

DESC....  02  FOR  CALCULATING  NG-C02  &  NG-EFF 

LL-LO  OP  LIMIT.  0.00 

LH-HI  OP  LIMIT.  20.00 

LR-RATE  LIMIT..  20.00 

NX-NEXT  BLK....  NG-C02 

CP-COPY  TO..... 


TAG . 0IL2-02 

DV-DEVICE . DUMMY 

HT-H/W  OPTIONS. 

lO-ADDR . 

SC-SIG  COND.... 

EL-LOEGU .  0.00 

EH-HI  BGU .  20.00 

ET-EGU  TAG . PCT 

OR-OUT  REVERSE.  N 
CS-COLD  START.. 

DESC....  02  FOR  CALCULATING  0IL2-C02  &  OIL2-EFF 

LL-LO  OP  LIMIT.  0.00 

LH-HI  OP  LIMIT.  20.00 

LR-RATE  UMIT..  20.00 

NX-NEXT  BLK....  OIL2-C02 

CP-COPY  TO . 


TAG . 0IL6-02 

DV-DEVICE. . DUMMY 

HT-H/W  OPTIONS. 

lO-ADDR. . 

SC-SIG  COND.... 

EL-LOEGU .  0.00 

EH-HIEGU .  20.00 

ET-EGU  TAG . PCT 

OR-OUT  REVERSE.  N 
CS-COLD  START.. 

DESC....  02  FOR  CALCULATING  0IL6-C02  &OIL6-EFF 

LL-LO  OP  LIMIT.  0.00 

LH-m  OP  LIMIT.  20.00 

LR-RATE  UMrr..  20.00 

NX-NEXT  BLK....  0IL6-C02 

CP-COPY  TO . 


TYPE:  TR 

TAG _ EFF-TR 

EL-LOEGU......  45.00 

EH-HIEGU .  95.00 
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ET-EGU  TAG . PCT 

AT-AVG/CX>MFRESS..  1 
NX-NEXT  BLK.... 


TYPE:CA 

TAG . OIL2-C02 

EL-LOEGU .  0.00 

EH-HIEGU .  100.00 

ET-EGU  TAG . PCT 

RN-ROUND  ENAB..  Y 
OUTPUT  =  (B-(C*A)) 

B- INPUT  2....  15.611 
C- INPUT  3....  0.7425 
D- INPUT  4....  0.0000 
E- INPUT  5....  0.00 
F- INPUT  6....  0.00 
G- INPUT  7....  0.00 
H- INPUT  8....  0.00 
NX-NEXT  BLK....  OIL2-EFF1 
CP-COPY  TO . 


TAG . Om2-EFFl 

EL-LOEGU .  0.00 

EH-HIEGU .  100.00 

ET-EGU  TAG.....  PCT 
RN-ROUND  ENAB..  Y 

OUTPUT  =  (((((B-(C*H)H(D*H)*H))-(E*A))+((F*A)*A))+<(G*A)*H)) 

B- INPUT  2....  202.08 

C- INPUT  3....  9.6511 

D- INPUT  4....  0.1732 

E- INPUT  5....  13.217 

F- INPUT  6....  0.4135 

G- INPUT  7....  0.5903 

H  -  INPUT  8....  02IUNK 

NX-NEXT  BLK....  OIL2-EPF 

CP-COPY  TO . 


TAG . OIL2-EFF 

EL-LOEGU......  0.00 

EH-HIEGU. .  100.00 

ET-EGU  TAG . PCT 

RN-ROUND  ENAB..  Y 
OUTPUT  =  (A-(B*Q) 
B- INPUT  2....  0.0268 
C  -  INPUT  3....  TEMP 
D- INPUT  4....  0.00 
E- INPUT  5....  0.00 
F- INPUT  6....  0.00 
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G- INPUT  7....  0.00 
H- INPUT  8....  0.00 
NX-NEXT  BLK.... 
CP-COPY  TO . 


TAG . NG-C02 

EL-LOEGU .  0.00 

EH-HIEGU .  100.00 

ET-EGU  TAG . PCT 

RN-ROUND  ENAB..  Y 
OUTPUT  =  (B-<C*A)) 

B- INPUT  2....  11.801 
C- INPUT  3....  0.5615 
D- INPUT  4....  0.0000 
E- INPUT  5....  0.00 
F- INPUT  6....  0.00 
G- INPUT  7....  0.00 
H- INPUT  8....  0.00 
NX-NEXT  BLK....  NG-EFFl 
CP-COPY  TO . 


TAG . NG-EFFl 

EL-LOEGU .  0.00 

EH-HIEGU .  100.00 

ET-EGU  TAG . PCT 

RN-ROUND  ENAB..  Y 

OUTPUT  =  ((((((-B)+(C*H)H(D*H)*H))+(E*A)H(F*A)*A))+((G*A)*H)) 

B- INPUT  2....  80.655 

C- INPUT  3....  16.757 

D- INPUT  4....  0.4514 

E- INPUT  5....  29.450 

F- INPUT  6....  1.2529 

G  -  INPUT  7....  -1.4427 

H  -  INPUT  8....  02JUNK 

NX-NEXT  BLK....  NG-EFF 

CP-COPY  TO . 


TAG _ NG-EFF 

EL-LOEGU .  0.00 

EH-HIEGU .  100.00 

ET-EGU  TAG . PCT 

RN-ROUND  ENAB..  Y 
OUTPUT  a  (A+(B*C)) 
B  -  INPUT  2....  -0.0277 
C  -  INPUT  3....  TEMP 
D- INPUT  4....  0.00 
E- INPUT  5....  0.00 
F- INPUT  6....  0.00 
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O- INPUT  7....  0.00 
H- INPUT  8....  aOO 
NX-NEXT  BLK.... 
CP-COPY  TO..... 


TAG . 0IL6-C02 

El^LOEOU 0.00 

EH-HIEOU 100.00 

BT-EOU  TAG . PCT 

RN-ROUND  ENAB..  Y 
OUTPUT  •  (B-(C*A)) 

B- INPUT  2....  16.495 
C- INPUT  3....  0.7848 
D- INPUT  4....  aOO 
E- INPUT  3....  0.00 
P- INPUT  6....  0.00 
G- INPUT  7....  0.00 
H- INPUT  8....  0.00 
NX-NEXT  BLK....  OIL6-EFFI 
CP-COPY  TO..... 


TAG . OIL6-EFF1 

EL-LOEGU .  0.00 

EH-HIEGU .  100.00 

ET-EGU  TAG.....  PCT 
RN-ROUND  ENAB..  Y 

OUTPUT  *  (((((B+(C*H))+((D*H)*H))+(E*A))+((F*A)*A))+((G*A)*H)) 

B- INPUT  2....  23.291 

C  -  INPUT  3....  -18.325 

D- INPUT  4....  1.0014 

E- INPUT  5....  14.700 

P  -  INPUT  6....  -0.6183 

G- INPUT  7....  0.8204 

H  -  INPUT  8....  02JUNK 

NX-NEXT  BLK....  OIL6-EFF 

CP-COPY  TO . 


TAG . OIL6-EFF 

EL-LOEGU......  0.00 

EH-HEGU .  100.00 

ET-EGU  TAG.....  PCT 
RN-ROUND  ENAB..  Y 
OUTPUT  « (A+(B*Q) 
B- INPUT  2.... -0.0272 
C  -  INPUT  3....  TEMP 
D- INPUT  4....  0.00 
E- INPUT  5....  0.00 
F- INPUT 6....  0.00 
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G- INPUT  7....  0.00 
H- INPUT  8....  0.00 
NX-NEXT  BLK.... 
CP-COPY  TO . 


TYPE:PG 

TAG . FUEL-MODE 

AO-ALM  DEST....  TH 

AP-ALM  PRI . L 

ST-SCAN  TIME...  1 

IS-INTT  SCAN...  ON 

lA-INTT  A/M....  AUTO 

DESC...  Choose  efficiency  calculation  models 

00  IF  FUEL  =  1  GOTO  6 

01  NUL 

02  IF  FUEL  =  2  GOTO  10 
03  NUL 

04  IF  FUEL  =  3  GOTO  14 
03  NUL 

06  SETOUT  NG-02  02JUNK 
07  SETOUT  C02  NG-C02 
08  SETOUT  EFF  NG-EFF 
09  GOTO  17 

10  SETOUT  0IL2-02  02JUNK 

11  SETOUT  C02  OIU-C02 

12  SETOUT  EFF  OIU-EFF 

13  GOTO  17 

14  SETOUT  0IL6-02  02nJNK 

15  SETOUT  C02  OE^C02 

16  SETOUT  EFF  OIL6-EFF 

17  GOTO  0 

18  NUL 

19  NUL 
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